u- 

rv'.orj 



■jL 



NPS-57GN71021A 



// 



United States 
Naval Postgraduate School 




DEVELOPMENT OF A UNIFIED HEURISTIC MODEL 
OF FLUID TURBULENCE 



by 



T. H. Gawain, D. Sc. and J. W. Pritchett 



February 1971^ 



Approved for public release; distribution unlimited. 



FEDDOCS 

D 208.14/2:NPS-57GN71021A 




NAVAL POSTGRADUATE SCHOOL 
Monterey , California 



Rear Admiral Robert W. McNitt, USN M. U. Clauser 

Superintendent Academic Dean 

ABSTRACT: 

A unified heuristic model of fluid turbulence is proposed v/hich supplements 
the basic equations of motion and continuity so as to define determinate 
solutions under a wide variety of conditions. A single set of equations 
applies to all cases. Only the boundary conditions change for each 
application. 

The Reynolds stresses are related to the mean flow heuristically , through 
an eddy viscosity. This also involves the turbulent energy, which is 
analyzed with the aid of empirical expressions for dissipation and diffusion. 

A characteristic local length scale is introduced, based on a generalization 
of von Karman’s mixing length. 

The model is illustrated by application to a channel and a jet. Agreement 
with experiment is satisfactory. 

The report also gives the detailed rationale which underlies the model, 
including an original analysis of the stress/s train rate relation. Suggestions 
are made for further research. 
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Fart I 

The Heuristic Model and Its Applicatio n 
by 

T. H. Gavvain, D.Sc. 
and 

J. V. Fritchct 
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1.1 TNTFODUCTIOr 



The Nava] Radiological Defense Laboratory, prior to its disestablish- 
ment in June 1969, had for some time beer pursuing studies in numerical 
hydrodynamics in the field of underwater explosions. It is v’ell l.nown 
that fluid turbulence forms an important aspect of this problem. In fact, 
the lack of an adequate theory of turbulence has been an important limita- 
tion on the development of numerical hydrodynamics for applications not 
only to the explosion problem, but also to fluid mechanics problems of 
every type. There continues to be an urgent need to overcome this obstacle. 

Consequently, an agreement was reached between NRDL and the Naval 
Postgraduate School some time ago for a joint attack on the problems of 
fluid turbulence. Two main approaches were investigated. The first approach 
emphasized fundamental understanding of the phenomena without undue concern 
about developing computational methods for application to practical problems. 
Ihe second approach emphasizes the development of a feasible, heuristic 
computational method at the expense, if necessary, of analytical rigor. 

This report documents the excellent results that have been achieved in 
connection with this latter goal. Reference (1), although nov; largely 
superseded by more recent work, documents some earlier effort along the 
other avenue of attack. Further work is continuing at the Postgraduate 
School . 

Unfortunately, the Naval Radiological Laboratory v:^as disestablished 
in the midst of this research. The abrupt cessation of NRDL activities 
made it virtually impossible to complete certain details of the 
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mathematical model in a satisfactory manner. It was therefore decided to 
hold up distribution of a final report until such time as these additional 
points could be worked out. Meanwhile, the press of other sponsored research 
compelled the investigators to lay this task aside for a considerable time. 

Fortunately, it was eventually possible to complete and revise the 
details mentioned above. Consequently, the present report, which describes 
the final model, is now being released for publication at last. Although 
the major part of the work described was accomplished some time ago, the 
model itself should continue to be of scientific interest for some time, 
owing to the basic and long term nature of the problem with which it deals. 

As an interesting historic note, it seems worthwhile to mention that 
the heuristic model developed in this report resembles in some respects 
the approach first suggested by Prandtl in 1945, Ref. (2). It also has 
some similarities to, and some differences from, the more recent work of 
Harlow and his associates. References (3) and (4). 

In incompressible flow, whether laminar or turbulent, the equations 
of momentum and continuity, along with the boundary conditions, theoretically 
suffice to establish completely the exact fluid motion. If the flow happens 
to be turbulent, however, the actual detailed motion becomes so complex 
that, although it is in principle determinate, its actual calculation 
would impose an overwhelming computational burden. 

The usual way around this difficulty is to average the equations of 
momentum and continuity. This results, of course, in an enormous 
simplification. Unfortunately, it also involves a significant and irre- 
trievalbe loss of essential information. Consequently, owing to the 
presence of the unknown Reynolds stresses which are created by this aver- 
aging process, the averaged equations of momentum and continuity do not 
in themselves comprise a determinate set. 
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In order to define a determinate solution, additional relations are 
required to fix the unknown Reynolds stresses. Unfortunately, these 
supplementary relations cannot be established from the original equations 
by any purely deductive process. For this purpose, supplementary empirical 
hypotheses are an unavoidable necessity. 

From another viev/point, it may be stated that the averaged equations 
of motion show’ the effect of the Reynolds stresses upon the mean flow. 
However, the reciprocal effect of the mean flow upon the Reynolds stresses 
is lost in the averaging process. Hence some adequate hypothesis must be 
found for representing this relation, at least approximately. 

For this purpose, a heuristic approach w’hich seems plausible is to 
postulate a relation between the Reynolds stresses and the mean flow which 
is analogous to the relation that is known to govern the viscous stresses. 
The analogue of the ordinary molecular kinematic viscosity is the so- 
called eddy kinematic viscosity. The problem becomes, therefore, to 
determine empirically the general law^ which governs this mean effective 
eddy viscosity at every space/time point in the flow field. 

The eddy viscosity presumably depends on a number of variables, one 
of the most important of which is the local kinetic energy of turbulence. 
Therefore, it becomes necessary to find the space/ time distribution of the 
turbulent energy. Fortunately, the governing energy equation can be 
deduced rigorously from the original equations of motion. However, the 
energy equation itself introduces two additional unknowns which can only 
be approximated in the same heuristic and empirical fashion as the eddy 
viscosity itself. The additioral unknovms are the rate of dissipation of 
turbulent energy into heat, and the rate of turbulent diffusion of energy. 

Theory and experiment both show that the eddy viscosity, and the 
dissipation and diffusion functions as well, depend not only on the 
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turbulent energy itself, but also on a local length scale parameter which 
can be associated with each space/ time point in the flow field. Von 
Karman was perhaps the first to point out how a physically meaningful 
characteristic length can be defined in terms of local space derivatives 
of the mean velocity at any point in the flow. In the present paper, the 
original approach of von Karman is further developed and refined. It now 
takes into account not only the velocity derivatives at the designated 
point itself, but also the values in the general vicinity of the point. 

By employing dimensional analysis, and by applying the available 
experimantal data, we finally obtain three empirical expressions which 
determine to a reasonable approximation the eddy viscosity, the heat 
dissipation and the turbulent diffusion, respectively. These expressions 
also involve the turbulent energy, the local length parameter, and the 
distance to the nearest fixed wall (if any). Of course, these empirical 
expressions are amenable to further investigation and development. 

In this way a single consistent and determinate set of eouations is 
established which applies in principle to any incompressible turbulent 
flow field. Only the boundary conditions differ for each specific 
applicati on . 

The application of this heuristic model is shovnn for two widely 
different cases, namely, a two-dimensional duct, and an axi-sym.metric jet. 
The distributions of energy and Peynclds stress are computed for each 
case and compared with experimental measurements. The agreement is 
satisfactory. 

These results, although limited in extent, strongly suggest that the 
proposed heuristic model of fluid turbulence is basically correct, at 
least in its main essentials. 
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This report is divided into two parts. Part I summarizes the final 
heuristic model and illustrates its application to two typical cases, 
namely, a channel flow and a free jet. Part II fives the corresponding 
detailed bacVground, theory and rationale. It tells vhy the final model 
is just what it is. It shows bow tlie various parameters of the model are 
related to fundamental aspects of turbulent correlation and spectra. An 
original analysis is presented of the general tensor relation betv^een 
stress and strain rate. This provides the basis for judging both the 
usefulness and the limitations of the eddy viscosity concept as employed 
in the heuristic model. The energy equation is derived from: first 
principles, then transformed so as to be practically useful. This part of 
the report also raises various questions of a fundamental nature which 
could be fruitful for future research. 
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1.2 SYMBOLS 



I.ocal dimensionless lenpth scale in jet 


Fef. Eos. 
1.6-31 


Overall dimensionless length scale of 
cross-sec tion in jet 

Cross-sectional area 


1.6-28 


Channel half width, or jet reference length 


1.6-1 


Functions in momentum equations for jet 


1.6-9 thru -14 



2 2 2 

Exponents in definitions of I , J , X , 
and T 


2.6-14 thru -19 


Arbitrary constants in definitions of 
A 2 and t 


2. 6-8, -9 


Determinant. Also diamieter of jet 


2.6-16 


Mean kinetic energy of turbulence per 
unit mass 


1.3-36, 2.3-5 


Instantaneous kinetic energy of turbulence 
per unit mass 


A- 12 


Rate of dissipation of turbulent energy 
per unit mass into heat. 


2.8-8 


Gaussian function 


3.6-6 


Ratio of local length scale of mean flov^ 
to correlation length 


2.5-1 


Ratio of local time scale of mean flow 
to correlation time 


2.5-2 



Non-isotropic tensor in p^eneralized stress/ 2.7-16 
strain rate lav^ 

Denotes values of f.. for principal axes 
of strain rate 



F 

F' 



Dimensionless stream function of jet 
Denotes derivative 



'’dF'' 



dn 



1 . 6-1 



F, G, H. 



S’ ^3 



,(S) 



H 



(T) 



i, j, k, 1 

2 2 
I , J 



Z 

(n) 



r, z, 0 



R. . 



R 



(S) 

ij 



F 



(T) 

ij 



Functions in energy equation for two- 
dimensional channel 

Parameter in stress/straln rate law 

Functions in analysis of jet 

Integral of weighting function 

Integral of spatial weighting function 

Integral of temporal weighting function 

Tensor indices 

Characteristic Integrals for a point 
Von Karman^s mixing length 
Superscript indicating realization number 
Reynolds number 

Kinematic pressure (pressure divided by 
density) of turbulent fluctuations 

Kinematic pressure (pressure divided by 
density) of mean flow 

Root mean square turbulent velocity 
fluctuation 

Cylindrical coordinates: radial, axial, 

and angular, respectively 

Correlation tensor 
Spatial correlation tensor 
Temporal correlation tensor 



1.5-25, -26-27 



2.7-28 



1.6-25, -26, -27 

1.3-11 

2.5-9 



2.5-10 



1.3-13 and -14 
1.5-7, -8 



1.3-16 



2.3-6 



2.3-9 



2.3-9 



s 



t 



Exponential constant in Gaussian velocity 1.6-6 
distribution for jet (s = 0.102) 

Time 
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At 

T 



'ij 



T. , 

t: . 

ij 



u . 
1 

U. 

1 

u: 

1 



U, V 



u 




V* 

dv’ 

W 



W 



(S) 






V 



w 



,(S) 



V7 



(T) 



— f 

X 

Ax 



Time increment 2.3-1 

Semi-interval for time averaging 2.2-A 

Viscous stress tensor for turbulent A.-16 

fluctuations 

Viscous stress tensor for mean flow 

Viscous stress tensor for mean flov; plus A-2 
turbulent fluctuation 

Turbulent velocity fluctuation 2.2-2 

Mean velocity 2.2-1 

Velocity of mean flow plus turbulent 
fluctuation 



Radial and axial components, respectively, 



of axi-symmetric mean flow 

Reference velocity of turbulent jet 1.6-1 

Kinematic Reynolds stress tensor (stress 2.7-1 

divided by density) 

Friction velocity 1.5-2 

Element of volume 

Weighting function 1.3-10 

Spatial weighting function 2.5-4 

Temporal weighting function 2.5.5 

Normalized weighting function 1.3-12 

Normalized spatial weighting function 2.5-12 

Normalized temporal weighting function 2.5-13 

Position vector of arbitrary fixed 
reference point 

Position vector of variable point 

Displacement vector 1.3-9 
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> 


Cartesian coordinates 


X, y, 




y 


Distance to nearest wall, or distance 
from centerline of jet or center plane 
of channel 


z 


Axial coordinate. Also nondimensional 

wall distance exponent 1.5-23 


a 


Dimensionless eddy viscosity coefficient 1.3-5, -18 


P, 6* 


Dimensionless heat dissipation coefficients 1.3-19,2.8-10 


Y 


Dimensionless diffusion coefficient 1.3-20,2.8-21 




Strain rate tensor of turbulent A-17 

fluctuations 




Strain rate tensor of mean flow 1.3-4 


r* r* r* 
11*^22’^23 


Strain rate tensor components for principal 
axes of strain rate 


6. . 

iJ 


The Kronecker delta tensor, with value unity 
for equal Indices, but zero otherwise 


e 


Effective eddy kinematic viscosity 2.7-22 


e 

o 


Eddy viscosity in centerline of jet 1.6-16 


^jk’^ijki 


Tensor versions of eddy viscosity 2.7-13,-14 


e" 


’‘Pressure diffusion" coefficient 2.8-17 


e’ 


Kinetic energy diffusion coefficient 2 #8-18 


; 


Radial coordinate B-H 


n 


Dimensionless distance from center plane 1.5-1 

of two-dimensional channel, or from 
centerline of jet 


CD 

<• 

CD 


Angular coordinate 


e* 


Strain rate parameter in stress/strain rate 2.7-26 
law 


K 


Constant in von Karman’s mixing length 1.5-8 

theory 
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K , K . 
1 



1 

A* 



D 



T . . 



t! . 

tV. 

IJ 



,2 



1.2 



ij 



0) 

Q 



Wave number vector 1.5-8 

Local length scale of mean flow 1.3-15 

Overall length scale 1.3-8 

Correlation distance 2.3-11 

Dissipation length 2.8-8 

Molecular kinematic viscosity 

Separation distance in meridional plane B-6 

Summation sign 

Local time scale of mean flow 2.6-18 

Correlation time 2.3-12 



Kinematic Reynolds stress tensor (stress 2.7-1 

divided by density) 

Distortional kinematic stress tensor 2.7-A 

Residual non-isotropic kinematic stress 2.7-19 

tensor 

Viscous distortional stress tensor 2.7-10 



Invariant of the tensor 

Invariant of the tensor tV. 

13 

Fourier transform of the correlation 

tensor F . . 

13 

Stream function for jet. Also auxiliary 
weighting function for axi-symmetric flows 

Stress parameter in stress/strain rate law 

Angular frequency or temporal wave number 

Invariant of strain rate tensor 
representing generalized shearing rate 

Absolute magnitude of the gradient of Q 



2.7-5 



2.7-20 

2.A-2 

1 . 6-1 

B- 13 



2.7- 29 



1.3- 6 

1.3- 7, 2. 



Note: In the jet analysis, subscript 1 on any quantity 1.6-33 

denotes the values at z = 1. It is therefore 
always a function of n only. 



1.3 SUmAFY OF BASIC EQUATIONS OF THF HEURISTIC MODEL 



The purpose of this section is to summarize the basic equations 
which, taken together, constitute the proposed heuristic model of fluid 
turbulence. A few explanatory comments arc offered to identify these 
various equations and to show their general relation to the overal] theory, 
but no detailed justification is given in this section; such details are 
included in other parts of this report. Also, the eouations in this 
section are sthted in three dimensional cartesian form. The reductions 
to forms specifically applicable to plane or axi-symmetrical mean flows 
are treated elsewhere. 

An incompressible mean flow field must everywhere satisfy the 
continuity condition 

^ = C (1.3-1) 

3x . 



and the ensemble averaged Navier Stokes equations of motion 





. 9 


f '' 

u.u. 
1 jj 




ra u. 1 

. . A . 


fip_l ^ f’-'j.ji 


^9t 


3x . 
J 


= V 


9x . 9x , j 
J 


[9xJ [ax. J 



(- 2 ) 



as well as the particular boundary conditions appropriate to the specific 
flow field under consideration. 

However, because of the occurence of the unknown Reynolds stress 
Eqs. (-1) and (-2) do not in themselves constitute a determinate set 
of equations. In order to define a determinate solution, additional relations 
are required. Purely deductive theory alone does not provide these relations. 
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Hence appropriate supplementary hypotheses must be tentatively postulated, 
subject to eventual verification, modification or disproof by experimental 
evidence . 

For the reason discussed at length elsewhere, we postulate that 
the Reynolds stresses can be adeauately related to the strain rates of the 
mean flow through a law which has been deliberately somewhat oversimplified 
to the form 



ij 



^ . + cr. . 

3 IJ IJ 



(-3) 



In this equation, q is the mean square velocity of the turbulent 



fluctuations, e is the so-called eddy kinematic viscosity, and T is the 
strain rate of the mean flow as defined by the expression 



r. . 

ij 



8x. 3x. 

^ 



(- 4 ) 



Unfortunately, the eddy viscosity c as defined above is not a simple 
property of the fluid itself, but rather a complex property of the turbulent 
flo\j field. Physical and dimensional reasoning leads to a supplemantary 
postulate, namely that c is expressible in the following form. 




(-5) 



Here X represents a local length scale, defined below, which varies 
from point to point in the field, but which is everyxr^here related in a 
consistent way to the local mean velocity distribution. The dimensionless 
coefficient a is a slowly varying universal function not predictable from 
theory, but determinable from experimental data. This function is specified 
later in Eq . (-18) 
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To establish the quantity A at any point x, we proceed as follows. 
We first compute the following mean flow variables throughout the field, 
namely. 



Q 



2 



1 

2 



r. .r, . 

ij 



(- 6 ) 



( 9 . 9 /)^ = 



'd9^' 




9x . 1 





(-7) 



Then we compute the reference quantity 



A 



2 

1 



/n^dv 

f(99l')^dv 



(- 8 ) 



where the volume integrals extend over the entire field. Of course is 
therefore a characteristic constant for the field as a whole and is 
independent of x, but it may vary as a function of time if the mean flow 
be unsteady. 

The calculation continues as follows. All integrals extend over 
the entire field of flow. Let 

X = coordinate of the fixed point at which X is being evaluated 

5c* = coordinate of variable point at wiiich integrands below are 
evaluated 

Ax= (x* - x) = separation variable (-9) 



W 



A X* Ax 




e 



= weighting function 



(- 10 ) 



K = /Wdv* 



w 




= normalizing integral 
= normalized weighting function 



(- 11 ) 

(-32) 
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(~ 13 ) 



2 4 

I = /wQ dv = first characteristic integral 
2 2 

J = fw{VcU') dv = second characteristic integral (~14) 

2 2 - 

Of course the integrals T and are functions of x and, if the 
mean flovv is unsteady, of time t as weJl. 

The required length scale A at any point x is finally given by the 

ratio 




The above calculation sequence establishes not only the desired 
local length scale A but also the useful auxiliary quantity J. It turns 
out that J is needed below in connection with the energy enuation. The 
rationale involved is fully explained elsewhere in this report. 

If the mean flow happens to be either plane or axi-symmetri c , all 
of the above volume integrals can be reduced to corresponding surface 
integrals. The reduction for the axi-symme trie case is rather interesting; 
it is given in Appendix B. 

If the flow field happens to be self-similar, the above calculations 
need be made on only one representative cross-section. Corresponding 
results for other cross-sections are then found simply by applying the 
proper scaling rules. If the flow field is not self-similar, all points 
must be computed independently. If the flow field is steady, the calcula- 
tions can be made once and for all. If it is unsteady, they must, of 
course, be repeated at appropriate intervals of time. 

Once the distribution of A has finally been established, it is then 
necessary to determine the distribution of the turbulent energy E, or 
equivalently, of the root mean sauare turbulent velocity 
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q = 



(-16) 



This is governed by the energy equation 






i i 




(>17) 



Of course the eddy viscosity e in the energy equation is the same 
as that earlier defined in Eq • (-5) in connection with the stress/strain 
rate law. If desired, Eq . (-5) may be used to eliminate e from the 
energy equation. In any case, the coefficient a which governs c is fixed 
by Eq . (~18) below. 

Each term in Eq . (-17) has dimensions of time rate of change of 
energy per unit mass. The four terms on the right represent, respectively, 
generation of turbulent energy, dissipation to heat, convection by the 
mean flow, and turbulent diffusion. Molecular diffusion is negligible. 

The dimensionless coefficients p and y which occur in the energy 
equation, like the dimensionless coefficient a introduced earlier, are 
assumed to be w^ell behaved universal functions. However, they are not 
determinable from theory, but must be established empirically. 

On the basis of the work done so far, the follov;ing expressions 
seem to provide satisfactory agreem.ent with the available data. 




(- 20 ) 



(-19) 



(-18) 
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In these expressions, y is the distance to the nearest fixed 



boundary. As the boundary is approached y approaches zero, but it turns 
out that under these conditions, the Jenpth parameter A also approaches 



arbitrarily be assigned a lower limit of unity. On the other hand, far 
from any fixed boundary, or in the absence of fixed boundaries, ve set 
y = «>, v;hereupon a, f, and y reduce to three simple constants. 

The proposed heuristic model of fluid turbulence consists, therefore, 
of the foregoing equations (-1) through (-20) inclusive, plus the particular 
boundary conditions appropriate to the specific problem under consideration. 
This system of equations theoretically defines a determinate solution for 
the general case of inhomogeneous and non-s tationary fluid turbulence. It 
fixes the mean flov/ field, the distribution of turbulent energy and the 
Reynolds stresses . 

The results obtained so far tend to support the proposed model, but 
further v/ork is necessary to ascertain more fully its adequacy and range of 
applicability. Various features of the model arc readily amenable to 
modification if necessary, so as to iuprove its range and accuracy. 

An important merit of the model, if it continues to prove successful , 
is that it reduces a v^ide range of dissimilar flow geometries and flow 
conditions to a single unified tb.cory. 



y 

zero simultaneously, so that the ratio ^ remains finite. In fact, it 



happens that at the wall itself, equals unity (not zero). In any case, 

it is stipulated in connection with the above formulas, that ^ shall 




1.4 CEbEPAI yiTHOb OP SObUTinN 



Although the equations (-1) through (-20) of the previous section, 



along with the appropriate boundary conditions, theoretically suffice to 
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dctcrnine all characteristics of the flovj field, the actual solution of 
this system of equations for any non-trlvial case is by no means simple. 

Let us consider first the general case of an unsteady mean flov^, with the 
initial velocity distribution of time t = 0 arbitrarily prescribed. 

Of course, the prescribed initial velocity distribution must satisfy 
the continuity condition, Eq . (-1). This offers no particular difficulty. 

In the usual case in which the mean flov; is either plane or axl -symmetrl c , 
the continuity condition may he satisfied identically by expressing the 
velocity components in terms of a stream function. Let us suppose that 
this has been done, and tliat the corresponding mean flow velocity components 
at time t=0 have been calculated. The sulseauent calculation procedure may 
now be described in somewhat general terms as follov/s. 

The next step is to perform the sequence of calculations summarized 
in Eqs. (-4) through (~15) . In this v/ay is obtained the local length 
scale > at every point in the field, and tlie quantity J. 

Next we must solve the energy equation (-17) for the mean souare 

turbulent velocity q at every field point. For this purpose v:e substitute 

Ka . (-5) into (-17) and use Eqs. (-18), (-19), and (-20) to evaluate the 

three empirical coefficients a, , and y. (In some cases the actual 

2 

distribution of the turbulent energy E = may be knoxv/n from experimental 
measurement. If so, v/e have at this point one important comparison between 
theory and experiment.) 

Once estimates of X and ci are available from the foregoing calculations, 
it becomes a simple matter to solve Eq . (-5) for the eddy viscosity e, and 
Eq . (-3) for the Peynolds stresses. 
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The Feynolds stresses obtained in this vray ir.ay now be substituted 
back into the bavier Stokes relations, Frjs. (~2) . By taking the divergence 
of hqs. (-2), \JC find that 



o 



3x . ox . 

1 3' 



k ^ [b.U. + T. J 

Bx . 9x . 31 31 

1 J 



(l.A-1) 



Since mean velocity and stress are at th.is point I nown, the nl ove result 
can be solved for the pressure distribution, subject to the specified 
boundary conditions on pressure. This part of the solution involves a 
lengthy iteration procedure, but is straigh tf orv^ard in principle. 

Once b., T.., and P are all established, we can revert to En . (-2) 

in . X / 



St ) 



.1 at 



every field 



and solve for the time rate of change of velocity 
point. By integrating forward tlirough time over a small interval, v.^e 
can find the new mean velocity distril. ution at the later time. 

The atove procedure may be repeated as often as desired, thus 
tracing the gradual evolution of the flow field over time. If the boundary 
conditions happen to be stationary in time, the entire flow field will 
evolve toward a stationary state, thus defining the solution of the 
corresponding steady state problem. 

The al'ove process of reaching a steady state asymptotically in 
time can be quite expensive in terms of computational effort and time 
required. It is often possible to shorten or abridge the above calculation 
in various x^^ays, depending on the situation. For example, instead of 
substituting the Feynolds stresses predicted by Fn . (-3) into the Navier 
Stokes lav; Fn . (-2), we can solve the latter for the Feynolds stresses 
"required'* in order to maintain the prescribed mean m.otion in a steady 
state. These "required" stresses can then he com^pared with the stresses 



'predicted" from Eq. (-3). If the agreement betv;een the two sets of 



19 



stresses is deemed satisfactory, this can be taken as a verification of the 
adequacy of the piven velocity distribution as a estimate of the steady 
state. This is the procedure that lias been followed in our sample problems, 
and it has proved entirely satisfactory for the intended purpose. In 
this case the purpose has been to demonstrate the degree of agreement of 
our heuristic mode] with experimental observations. The agreement is 
considered satisfactory for both cases analyzed. 



1.5 APPLICATION TO A TWO-DIMENSIONAL CHANNEL 



Consider the steady, incompressible, turbulent flow through a two- 
dimensional channel of uniform height 2b as shown in Fig. 1.5-1. Let y 
be distance from the channel center-line. Choosing b as the reference unit 
of length, the corresponding dimensionless distance becomes 




Let the kinematic shear stress at the wall be 



(1.5-1) 



(- 2 ) 

where v* is the so-called friction velocity. Without loss of generality, 
we can always choose reference units of time such that 




The kinematic shear stress x varies linearly across the channel so 
that at any location n we have 

T = n (-A) 

According to our present notation, if U is the mean velocity at any 

n , then 
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y,T? 




FIGURE 1.5-1 FLOW IN A TWO-DIMENSIONAL CHANNEL 
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= U' = 



(-5) 







= U'' = Q' 



According to the von Karman mixing length theory 
2 2 

T = Z Q 

where the mixing length i is given, in our notation, by 

2 2 
Si = K 






,2 



(- 6 ) 



(-7) 



(- 8 ) 



and where k is a constant which remains to be determined. 

Upon combining Eqs. (-A) , (~7) , and (-8), taking the square root 
of both sides ard rearranging, we obtain the basic relation 



P' 



P 



2 




(-9) 



This equation is readily integrable and gives the well known von 
Karman velocity profile vzhich in the present notation becomes 



\ = l; = U + - {in(l - /n) + /n) (-10) 

V’^ O K 

It is known that this results gives satisfactory agreement with experiment 
if we place 

K = 0.36 (-11) 

The resulting quantities of interest in the present analysis are 



= i- = i 

“ *• 2k(1 - /T?) 

and 

fi' = U" r 

n(l - /Ti) 



(- 12 ) 



(-13) 



22 



We next require the "averaging distance" defined by the equation 



>2 , . ^ 0 
X^ = iim 



^ ri 



Q dn 



Jo 



(- 14 ) 



(on')^dii 



Upon substituting expressions (-12) and (-13) into (-14), we find 
that in the limit as approaches unity, both integrals approach infinity, 
hut the above ratio of these integrals approaches zero. Hence v^7e obtain 
the unusual result that, in this case. 



= 0 



(-15) 



Now the local length scale X may be computed from the relation 



1 Ax*Ax 



e X 



X = £im 

X^H.0 



o 






1 A X* A X 



e A 






(-16) 



This result is rem^arkable! It shov;s that our characteristic 
length X is in this case proportional to the classical von Kaman mixing 
length i. Form Eqs. (-8), (-12), (-13) and (-16) ve find that 



^ 2 . 9 /-N 

A = - = -^ = 2/n(l - /n) 



(-17) 



Similarly, because of the exceptional result (-15) , the integral 



for J reduces simply to 

ri 



J = Him 
A ^-^0 



A X. A X 



e ^ (Ofi’)^dn' 



1 _ Ax* Ax 
0 

^1 

e dn' 



= (nn')^ 



(- 18 ) 
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Therefore again using (-12) and (-13) we find that 



j = nn' = 



8K^/n(l - /n)^ 



(- 19 ) 



The three empirical functions a, and y may now be expressed in 
the form 



= a = 0.065(1 + e ^ } 



A JlV. 



(- 20 ) 



and 



where 



= 3.7(1 + e ^ } 



= 1.4 - 0.4 e 



-z 



2 if 1 ,12 

" ■ 4[tS - M 

With the functions 9., X, J, a, 6, and y now knovm, we can proceed 
to solve the energy equation. The unsteady term and the convective tenri 
both vanish from this equation in the present instance and we obtain 



(- 21 ) 



(- 22 ) 



(-23) 






Dn 



ayXv(2E 



9n 



= 0 



(-24) 



Energy E is now the only unknown. To simplify the equation further, 
we introduce the following notation. 



2 


= F(n) 


(-25) 




= G(n) 


(-26) 


ayX 


= H(n) 


(-2.7) 



Of course, these three functions are now known. Therefore the energy 
equation reduces to 
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(- 28 ) 



This equation cannot be solved analytically, but it can easily be 



expanded and rearranged for numerical solution by finite differences. 



The desired boundary conditions are 



At n = 0 [ 4 ^] = E' = 0 



(-29) 



At n = 1 E 



1 



0 



(-30) 



This last condition cannot be satisfied directly, however, we can integrate 

from the centerline out, using a trial value of the centerline energy E^. 

By trial and error, a value of F can then be found which will in fact 

o 

yield the required condition E^ = 0 at the wall. 

This is the method v?hich was used to compute the energy distribution 
shown by the solid line, in Fig. 1.5-2. It may be seen that the curve 
obtained in this way is in reasonably good agreement with the experimental 
data of Peichart and of Laufer, Refs. 5 and 6. 

Once the energy distribution is known, the dimensionless Reynolds 
stresses may be computed from the relation 



Except very near the walls, the above Reynolds stress should agree 
V7ith the total shear stress defined by Eq . (-4). 

If we wish to switch back to absolute units, Eqs. (-4) and (-31) 
reduce respectively to 



-uv 



T = aXs/2EQ 



(-31) 



1 = X 



Total Stress 



(-32) 




Reynolds Stress 



(-33) 
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Turbulent Energy E/v*^ 
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Measurements 
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Heuristic Model 
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X 
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FIGURE 1.5-2 TURBULENT ENERGY DISTRIBUTION IN A 

TWO DIMENSIONAL CHANNEL 
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Shear Stress -uv/v 




FIGURE 1.5-3 REYNOLDS SHEAR STRESS DISTRIBUTION 
IN A TWO DIMENSIONAL CHANNEL 
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In Fip, . 1.5-3, the solid line shows the above heuristic estimate of 
the Reynolds stress, computed for a Reynolds number of 10,000. The dashed 
straight line represents the true total stress, which should agree closely 
with the Reynolds stress everywhere except in the immediate vicinity of 
the wall. The agreement is seen to be good; the computed results are on 
the whole more accurate than are the experimental measurements from Reference 
(6), as shown on this figure. 

We conclude that this example tends to substantiate the proposed 
heuristic model of turbulence. 

1.6 APPLICATION TO AN AXI-SYMMETRIC JET 



Consider a free turbulent jet discharging into a quiescent atmosphere 
as shown in Fig. 1.6-1. The radial and axial coordinates are r and z and 
the corresponding velocity components are U and V, respectively. 

Experimental results indicate that the velocity profiles of the 
mean flow at various cross-sections are self-similar. Specifically, if 
z be measured from, a suitable virtual origin, the stream function ip can 
be reduced to the form 

= U bzF(n) (1.6-1) 

o 

where 




and where U and b are constants of the jet. 
o 

Consequently the velocity components of the mean flow become 
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FIGURE 1.6-1 THE AXI-SYMMETRIC JET 
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(- 4 ) 




where 



F’ = 






(-5) 



The known experimental results also indicate that the generalized 
velocity distribution through the jet Is very well approximated by the 
simple expression 



II 

n 



f 




(- 6 ) 



Where s is a characteristic constant for a turbulent jet* It has the known 
value s = 0.102, 

Furthermore, it is useful and convenient to adopt units of length 
and time such that, wh.en expressed in these units, * 1 and b = 1. 

Of course, tdiere is no loss in generality involved in this. It amounts 
merely to non-dimensionalizing all quantities with respect to and b as 
reference parameters. Hence symbols like e, p, u, v and so on will now 
represent the corresponding dimensionless versions of kinematic viscosity, 
kinematic pressure, velocity components and the like. 

The next step in the analysis is to combine the momentum equations 
(-2) and the stress law (~3) of section 1.3. The results may he transformed 
to the z, n coordinate system defined above. The foregoing dimensionless 
functions e and p may be introduced therein, and the results reduced. 

This reduction, although very lengthy, is s traigh tf orvjard . In this \vay 
the following results are obtained. 

np’ -h 2p-h A^(n)e* + A^Cn)^ = A^(n) (-7) 

"P* 4* 0 + B^(n)e* + B 2 (n)e = £ 3 ( 0 ) (- 8 ) 



30 



The and B’s are known functions defined as follows, where f is the 

Gaussian function defined by Eq. (-6). 



A^(n) 



_ 2 _ 

2 

s^ 



2n 

1 



J 



r f 



( 










fs - -2-1 


2n^ 




LI 




s2| 


2 

s 


4 

s J 



/jCn) 



(n) 



= f(l - 2f) 



1 _ 

s2 



r-f + --2 a - f) 
m" 



B^(n) = 



= 




- (3 - i,- 




V s J 





4ri^ 

72 “ J 



nf 



B3(n) 



= -qf + 



s 

n -f --- 
2n 



f(i - f) - a - f)2 

4 ri 



Tlieoretically , it is possible to eliminate p and p' between the two 
nomenturr. equations (~7) and (-8). This yields a single vorticity ecuation 
of second order in e v/hich can Le solved for ti e uni ro\m eddy viscosity 
distribution e (n) • 

liowever, it is known from exj eriment that the terms in p and p’ are 
negligible in Eqs. (-7) and (-8). I'oreover, the terms in In. (-8) are 
negligible in comiparison with tb.ose in i:q . (-7). Vie can tal e advantage of 
this information to simplify the solution for e(n). We simply discard 
Eq . (-8), and drop the terms in p and p* from Lq . (-7). The resulting 
equation is now only of first order in e . It may be v^ritten in the form. 



! 

/ 


(' ■ V 

^ S ' 


- 


[5 _ 

r s2J 


1 2n , 4 n 

1 2 4 

s s 


. . _ (1 - 2f) 




fl 2 . 2n2] 

1 - 2 + 2 

^ s s ^ J 


r b ' 


0 -,2'! 

, 2 2n 

i 2 2 

'■ S S'* 



- 9 ) 

- 10 ) 

-11) 

- 12 ) 

-13) 

*1A) 



- 15 ) 
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by 



Solution of this equation analytically is not feasible, but solution 
numerical integration is simple. 

bote first that at the centerline r, = 0, this ecuation pives at once 



e 

o 




0.00260 



(- 16 ) 



This result provides the boundary value required for the numerical integration. 

It is also of interest to note that for very large values of n , 

Eq . (-15) simplifies to the following approximation, namely, 




The solution is 




(-17) 



e 




-5 -sn , 

n e dn 



(-18) 



The true numerical solution approaches the above analytical curve asymptoti- 
cally at large values of p. 

VJith the dimensionless eddy viscosity distribution c(n) Vnown, it 

becomes a simple matter to find the corresponding dimensionless Reynolds 

shear stress t from the relation 
rz 





[ 3U 




'av' 




T = -UV = C 










rz 




1 







By making use of Eqs. (-3) through (-6) this result can be reduced to the 
following form for computation 



(-19) 




The shear stress computed from this eouation is shown by the dashed line 
in Fig (1.6-2) . 



(- 20 ) 
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To form the F.ean flow parameters required for the heuristic model 



v;e note that, in cylindrical coordinates, 



^ f 



+ 2 



f \ 2 



3z 









and 








2 


f 


2^ 


[k, 




[Sz ^ 





The change to the new coordinates ri, z is quite lengthy, but not 
difficult otherwise. The results finally obtained can be summarized in 
the following form. 

2 1 

r/ = G^(n) 

z 



(ra’)^ = ^ 



10 



GjCn) + (' 






vdiere the G*s are kno\^ functions defined as follo\7S. 



^^(n) = 1^; 



12 + 1 ^ + ~ - 15 

s s 



' 4 [2 




4 , 4n^ 


2 2 I 2 

^ S S ^S 


- 3 


n + 

s J 



+ 21^ - 2 

‘•ri 



f(l - f) + -^-^(l - f) 



G2<n) 



^ - 39 
2 2 
■ s s 



_ 2 _ . 12 _ 1 _ , , 

2 [ 6 4 2 
s '■ s s s 



- 3 



5 16 

" "“6 

s 






s 



- 2f) f 



(1 - 2f)f + ^(1 - f) 
n 



(- 21 ) 

(- 22 ) 

(-23) 

(-24) 

(-25) 

(-26) 
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G (n) 
3 



= [• 



- -J 



12+36 

4 2 

s s 



2 2 

r, + — 



4. 12 4. 1 oaI ^ 

— - + — 7 + — ;; - 2m j n 



s s 



f4 



'•s 



- 1 



6.8 8 

n + “z n 



+ 2(1 - 



.2 



(1 - 2f) f 



(- 27 ) 



- 4 






'■n 



2 r 1 

(1 - f)f - Vf2 - ' 



(1 - f)‘ 



Kote that these three functions are all finite on the axis n = ^ 
and that they all vanish at infinity. 

To find the avera^^ing distance in a non-dimensional form appropriate 



to this problem, we evaluate the expression 



/ OO 4 (»00 2 

rdr / G, (n)ndn 

^ ' 0 1 

’ r^iVcV.')^Tdr /" [G^(n) + G^(n)]r,dn 



(-28) 



It should be observed that both integrals in this last equation are 
evaluated from n = 0 to n = " over the surface z = constant, and that 
in this case the quantity z cancels identically from the right side of 
the equation. There is no integration with respect to z. The dimensionless 
averaging distance a^^ is constant and independent of z. 

The integrals involved in Eq (-28) can easily be evaluated numerically, 
thus fixing a^. It is found that a^ is roughly the same as s, the jet 
width parameter; this seems reasonable. 

2 

A similar procedure is involved In determining the Integrals I and 

2 

J and the resulting local length scale parameter. Also the effects of 
polar symmetry must be taken into account in the manner of Appendix B. 

Again the integrations are from n = 0 to n = ® over the surface z = 
constant, and again the variable cancels from the final result. Thus 



3A 



I = 



.OO 2 . 

1 


[n'-n]2 

1 a, J 

yiT’ i|- 


[•/nTT'l 


dr' 


-1 


fn '-ri 
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/.• 4^1 

(. ^2 j 


dn’ 

i 



8 



(-29) 



=-. 



:ehi^ 



10 



/: 



I ^1 



Jj(n) 

10 



/n’ 



Hence the dimensionless loca] lenrth scale becomes 



O 1 

= a^(n) = — ; 



ri2 



'•j' 



lj(n) 

J^(n) 



(-30) 



(-31 



2 2 

The integrals I^(n) and must of course le evaluated numerically. 

It turns out that the local length scale a(n) obtained in this may is a 
slowly varying function. Tlie value remains nearly constant at about 0.1 
across most of the jet, and declines only slo\'ly at relatively large 
distances from the centerline. 

The energy equation for the axi-syn>me trie turl ulent jet may be 
vriLten, first in cylindrical coordinates, 



oh 

at 



= 0 = el 



The various functions which appear in the energy enuation can nov’ 
be represented as functions, of the similarity coordinates z and r. as 
follows 



(“'12) 
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f = c(n) 



L = 



V = 



V = 



Ej^(n) 



U^(r,) 



V^(n) 



-(r) 



J^(n) 



(- 33 ) 



= zn 



l~( ) = ) 

dr z 0 r) 



Upor substituting the expressions(-33) nnto the energy equation, v:e 
obtain , 



^ _ g i 1 



1 ^ r TT p 1 

4 d n 11 d Z 



z n 



V F 
1 1 



L Z 



(-34) 



^ ^L. 

4 3n 
z n 



fd}\. 



r\ye 



9r> 



9z 



Yc 



rF ^ ■'» 

9 f 1 



9zl 2 
z 



= 0 



r r 



This result confirms tlie fact that all terms of the energy equation 
-4 

are proportional to z and that the above energy lax; is indeed compatible 
with the self-similarity of the overall solution. In order to expand 
this so as to be able to eliminate z from the result, xv"e male use of the 
relation 



d 

9z 



[ 1 



(-35) 



In this way we obtain the general result 



cpf - + 3VjEj + „ |;^IVjEjl 



n 9n 



f’hY 
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^ ^ — 1 , 



(-36) 
= 0 
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Considerable simplification can be achieved with little loss of 



accuracy by neglecting the effect of the slight divergence of the jet 
on the energy equation. This amounts to the assumption that 



iiJ 



0 



(- 37 ) 



In this case the energy equation simplifies to 



-2 ,7/6. 1/3 1 9 r „ 1 3 

er> - 6(2r. ) J ~ “ Tit'll’ ] + T TX 

1 1 1 HwH 11 pen 



— 




r\ye 


[anIJ 



= 0 



(-38) 



To complete the solution it is necessary to huow the empiricai 
functions a, 6, and y as defined by Eqs. (-18), (-19), and (~20) of 
Section 1.3. In this case, however, there arc no fixed walls involved, 
and these quantities become constants, namely, 
a = 0.065 

^ " 'TJ 

y = 1.4 



(-39) 



The eddy viscosity c(n) which occurs in the energy equation may be 
found in either of two ways. Firstly, we may use the solution previously 
found from the momentum analysis, by numerical integration of Fq . (-15). 
Secondly, we may use the basic postulate of the heuristic model, that is. 



e = aX/2E = aa72L^ 



(-40) 



This second alternative represents a somewhat more severe test of 
the heuristic model and was adopted for the numerical example here 
discussed. Thus substituting Eq . (-40) into (-38) gives 



aa/2E^r:!^ - P(2L^) 



1 ^ r IT r T ^ ^ ^ 
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(-41) 
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The energy function is the only unknoxvm remaining in this 

equation. The boundary conditions are 



At 


n = 0 


/9E >1 
^ = 0 
l,3ri j 


(-42) 


At 


n-K» 




(-43) 



Eq . (-41) must be solved numerically, 
for this purpose to use the substitution 




It happens to be convenient 



(~44) 



The equation is then expanded and the derivatives replaced by appropriate 
finite difference approximations in a straightforward manner. 

The boundary condition (-43) at infinity must be satisfied indirectly. 



Actually, the integration starts at the centerline using 



3q 



= 0 and 



a trial value of q^. A value of q^ is found by iteration which yields 
q ^ 0 at some large value of q. In this case the value n = 0.25 may be 
regarded as ‘'large.'* It is found that the precise position of the outer 
boundary makes virtually no difference in the results as long as it is 
greater than about 0.20. 

The energy distribution found in this way is shown by tlie solid line 
in Fig. 1.6-2. The agreement with the experimental results of Corrsin, 

Kef. (7), Laurence, Ref. (8), is seen to be satisfactory, especially in 
view of the degree of scatter in the data points. 

Once the energy distribution is knov/n the eddy viscosity is fixed 
by Eq . (-40) and the corresponding shear stresses by Eq . (-20). These are 
shown by the solid line in Fig. 1.6-3. The dashed line shows the corres- 
ponding shear stresses as computed from the momentum equation. Data points 
are from Corrsin, Ref. (9). 
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Turbulent Energy (E/Uo^XZ/b)^ 




Relative Radius r/z 



FIGURE 1.6-2 TURBULENT ENERGY DISTRIBUTION IN 
A CIRCULAR JET 
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Shear Stress (-uv/Uo^Xz/b)^ 




Relative Radius r/z 

FIGURE 1.6-3 REYNOLDS SHEAR STRESS DISTRIBUTION 
IN A CIRCULAR JET 
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If the heuristic model were completely correct, and if the assumed 
Gaussian velocity profile were exact, then the two shear stress curves 
shown in Fig. 1.6-2 would coincide exactly. As it is, the degree of 
agreement attained is considered to be reasonably satisfactory. The 
discrepancy between the two curves is smaller than the discrepancy between 
the experimental points and either of the curves. 

It seems probable that in this instance the theoretically computed 
values are actually more accurate than the experimentally measured ones; 
the experimantal measurement is innately difficult and uncertain. 

It may be concluded that, on a whole, these results for the turbulent 
jet substantiate quite well the proposed heuristic model of fluid turbulence. 



1.7 CONCLUSIONS AND RECO?I>^ENDATIONS 



It is concluded, on the basis of the evidence available so far, that 
the proposed heuristic model is basically adeouate for determining the 
principle flow characteristics in the general case of inhomogeneous and 
non-stationary turbulence in incompressible flow. 

It is recommended that the present model be applied also to other 
cases including pipe flow, two and three dimensional wakes, boundary layers 
and to some examples of unsteady flow. The aim would be to refine the 
present unified theory and extend its range of applicability. 

It is also recommended that further experimental information be 
obtained in connection with those aspects of the model for which the 
present data are insufficient. These aspects include, for example. 
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(a) the relation between the local length scale X of the mean flow 
and the correlation length X* of the turbulence, 

(b) the generalized three dimensional stress/strain rate relations 
which actually exist in regions of strong anisotropy such as 
in the flow near a wall, and 

(c) the influence of various key parameters on the rate of dissipation 
of turbulent energy into heat. 
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Part II 



Theory and Background of the Heuristic Model 

by 

T. H. Gawain, D.Sc. 
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2*1 THE FUNDAMENTAL PROBLEM 



Consider a general turbulent flow field, one which is both in- 
homogeneous and non-stationary . How can the turbulent structure of 
such a field be adequately analyzed and described? 

In principle, it is possible to formulate a valid general method 
for computing numerically the detailed solution u^^^(x,t) for a particular 
realization of the field, subject to specified boundary conditions and 
initial conditions* Once these conditions are fully specified, the 
equations of continuity and motion suffice to establish a fully determinate 
and detailed solution. In theory, this solution procedure could be repeated 
very many times to generate a whole ensemble of such detailed realizations. 
All realizations would be compatible with the same general macroscopic 
boundary conditions and initial conditions. However, they would differ 
more or less randomly in regard to the fine scale details of the boundary 
and initial conditions, and therefore also in the fine scale details of the 
resulting motions. 

Once such an ensemble of detailed solutions were available, any 
desired statistical features of the turbulence could then be computed 
from this ensemble in a completely determinate way. 

It is significant that the foregoing method does, in principle, 
satisfy the mathematical requirements for a determinate solution. This 
solution is based solely and strictly on the general equations of contin- 
uity and motion (the Navier Stokes equations), and on the pertinent bound- 
ary and initial conditions for the particular case in question. 
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In practice, however, the volume of detailed calculations required 
to carry out the above hypothetical scheme of calculation for any case 
of practical interest is simply overwhelming. It lies vastly beyond the 
capabilities of any computer yet known or contemplated. The basic reason 
for this, of course, is that the detailed space/time structure of turbulence 
is extremely fine grained and enormously intricate. Hence for all practical 
purposes, the above hypothetical solution procedure must be ruled out. 

There is a classical m.ethod for attempting to circumvent the above 
difficulty. Since it is impractical to solve the detailed equations and 
average the results, the classical method averages the basic equations 
themselves. It then attempts to compute the desired average characteristics 
of the flow from these averaged equations. Now, if the basic equations 
were linear, this technique v;ould indeed solve the problem. Unfortunately, 
they are not. As a result of the non-linearity, the averaging process 
introduces additional unknovms , the Reynolds stresses, into the basic 
equations. Consequently, there are finally more unknovms than available 
equations. Therefore, the averaged equations do not comprise a closed 
and determinate set. This constitutes the so-called closure problem of 
turbulence theory . 

It should be emphasized that the deficiency which arises from 
averaging the equations is f undamiental . No amount of mere manipulation 
of tlie equations, however involved, can remedy this defect. The averaging 
of the equations leads to a definite loss of essential information in 
comparison with the information implicit in the original equations. 

Therefore, so long as v:e are constrained to work only with the 
averaged equations instead of with the originals, it becomes unavoidably 
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necessary to invoke additional hypotheses in order to compensate for the 
lost information and to reestablish a determinate set of equations. 
Unfortunately, there is no known way to establish the required hypotheses 
deductively from the original fundamental equations themselves. Consequently, 
it is initially far from clear just what auxiliary hypotheses are most 
nearly adequate and satisfactory. 

Obviously, the basic justification for any particular hypotheses 
finally adopted must therefore rest on a comparison of the consequences 
deduced from these hypotheses with experimental observations. This is 
the accepted scientific method of inductive inference. While it is to be 
regretted that the problem apparently cannot be solved on a wholly 
deductive basis, there are no a-priori reasons for despairing over the 
prospects for ultimate success by the inductive method. 

It should also be pointed out at this juncture that the averaged 
equations of motion describe exactly the effect of the Reynolds stresses 
on the mean motion. However, they give no information whatever concerning 
the reciprocal effect of the mean motion on the Reynolds stresses. Yet 
some definite reciprocal effect of this kind unquestionably does exist. 

This, therefore, is precisely the problem with which our auxiliary 
hypotheses must necessarily deal. 



2.2 STATISTICAL DESCRIPTION OF TURBULENCE 



Consider an arbitrary and general flow field which is neither 
homogeneous in space nor stationary in time in regard to its turbulent 
structure. Let u^^\x,t) denote any one particular realization of this 
field. This means that the velocity components u^'^^'x, t) represent one 
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particular detailed solution of the equations of continuity and motion; 
this solution is compatible with the specified macroscopic boundary 
conditions and initial conditions. There exist, of course an infinity of 
such possible individual realizations which satisfy these general 
conditions. The separate realizations differ from one another only with 
respect to their fine scale details. We denote each realization by a 
distinct value of the index n, which is taken to be an integer. Thus we 
may have, for example, n = 1, 2, 3,...N where N is some indefinitely 
large number. Such a very large and representative collection of individual 
realizations we term an ensemble. 

Fortunately, the particular details of any one individual realization 
are seldom of direct interest. What is usually required is statistical 
information concerning the principal features of the ensemble as a whole. 

In particular, with respect to any variable of interest, what we would 
generally like to know are not the individual values for particular 
realizations, but rather the following items of information. These are 
listed in order of decreasing urgency. Thus 

(a) The average value of the variable over the ensemble. 

(b) The statistical variance of the quantity with respect to the 
ensemble average. 

(c) The detailed probability distribution of the quantity over 
the ensemble. 

( 

For example, if is the instantaneous total velocity at 

space/time point x,t for realization n, we require the corresponding 
ensemble mean 

U (x,t) = liim [ U!’'”'’0?,t) (2.2-1 

N-«o n=l 
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This result defines the mean flow field. Of course, the turbulent 
fluctuations of the nth realization are defined as the deviations from 
this mean. Thus 

u^''hx,t) = U|^"hx,t) - U^(x,t) (-2) 

From this definition, it follows that the turbulent fluctuation bas 

zero mean over the ensemble. 

The corresponding variance of this velocity fluctuation is given by 

-5 1 ^ 

u^(x,t) = £im I Oc,t) (No sum on index i) (-3) 

Finally, the probability function which defines the detailed statistical 
disttributlon of this random variable u^ over the ensemble may be denoted 
by the functional (x , t) ] . Here the quantity H^du^ represents the 

probability that for a realization n, drawn at random from the ensemble, 
the velocity component u^ at point x,t will have a value which lies 
between u^ and (u^ + du^) . If the space/time coordinates x,t be held 
constant at prescribed values, 4^(u^) reduces to a corresponding ordinary 
probability density distribution. Hence the functional ^ amounts to an 
ensemble of probability distributions, one for each space/time point of 
the field. 

If the turbulence happens to be stationary with respect to time in 
its statistical features, then it becomes permissible and convenient to 
replace the ensemble average by a simple time average over a single 
realization. The realization index n may therefore be dropped. For 
example, the mean velocity distribution may now be redifined as follows. 

+T 

U (S) = ilim ^ / U'(x,t)dt (-4) 

^ T-Mo -T ^ 
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Similar changes apply to the definitions of variance and probability 
distribution. 

Likewise, if the turbulence is homogeneous with respect to some 
line, surface or volume, the ensemble average may be replaced by the 
corresponding spatial average. 

In the following analysis, the average will always be an ensemble 
average, unless specified otherwise. Also the discussion will center 
largely on certain important averages and variances, but will not deal 
specifically with the corresponding detailed probability distributions. 



2.3 CORRELATION FUNCTIONS 



In order to formulate an adequate theory of turbulence, it is 
necessary to define certain basic concepts which characterize the 
turbulence in the general vicinity of any prescribed space/time point. 

In this connection, the two point velocity correlation tensor is of 
basic importance. 

Let P denote a reference space/time point with coordinates Ic , t 
regarded as arbitrary but fixed. Let P* and P'’ be two associated points 
with coordinates 

X* = X + Ax 
t* = t + At 

and (2.3*1) 

x” = X * Ax 
t” = t * At 

The quantities 2Ax, 2At represent the space and time separation between 
the points P’ and P” and are regarded as variable. Note, however that 
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P' and P” are always disposed symmetrically equal and opposite with 
respect to the fixed reference point P. 

Components of the turbulent velocity fluctuations at the variable 
points P* and P*' are denoted by 



and 



u^ = u^(x 4- Ax, t + At) at P* 



u” = ^ 



The mean product of these two velocity components is the tensor 



(- 2 ) 



u’u'! = u.(x + Ax, t + At) u. (x - Ax, t - At) 



i j ' ' J 

where the overbar denotes an ensemble average. 

At zero separation, the above tensor reduces (upon reversal of 
algebraic sign) to the ordinary Peynolds kinematic stress at the reference 
point P. Thus 



(- 3 ) 



T.. = -U.U, = -U, (X,t)u , (x, t) 
ij 1 J k ’ j ’ 



(- 4 ) 



By reason of its symmetry the stress tensor most 



six distinct components. The general velocity correlation tensor 
however, is not necessarily symmetrical with respect to an interchange 
of the indices and may therefore have up to nine distinct components. 

Upon contracting the Reynolds stress tensor, we obtain the 
physically important quantity 

q^(x,t) - u^u^ ” '^l ^ ^2 ^ ” 2E(x,t) 



(- 5 ) 



here E represents the ensemble mean turbulent energy at the point x,t and 
2 

q is the corresponding mean square velocity. 

We now define the space/time correlation tensor merely by 



normalizing with respect to the value of q at the fixed reference 

point P. Thus 
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i 



P. . 




(x, t;Ax,At) 



u.u. u , vx -r £ix, t + At)u.(x - Ax, t - At) 

1 3 ^ _i ! J ^ 

2 2 
q q (x,t) 



(- 6 ) 



It is of interest to compare the above definition which we have 
adopted with certain alternative forms which are sometimes encountered. 
One such alternative is 



u.(x,t)u.(x + Ax, t -b At) 

P! . = — ^ 

q(5T,t)q(x + Ax, t + At) 

Another alternative is 



(-7) 






u^(x,t)u^(x -t" Ax, t -K At) 

^ u^(x,t) (x ^f^Ax, t + At) 



(No sum on i or j) 



(- 8 ) 



The form we have chosen to preferable because it is related in a 
simpler fashion to the energy spectrum of turbulence, and because it has 
certain useful symmetry properties which will be pointed out presently. 

Two important special cases of the foregoing are the spatial and 
temporal correlation tensors, defined respectively as follows: 



/V u . (x + Ax , t)u . (x - Ax , t) 

- R (3c,t;Ax,0) - ^ ^ 

q (x,t) 



(-9) 



u (x,t + At)u (x,t - At) 

= R,.(x,t;0,Ax) = r J 

^ ^ q (x,t) 



Of particular significance are the first invariants of these two 
tensors, defined as follows: 
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il il 



R (x,t;Ax,0) - — 



u.(X + AX, t)u (X - AX, t) 



q^(x,t) 



(- 10 ) 



“ R^^(x,t;0,At) 



u^(x,t + At)u^(x,t - At) 



q^(x,t) 



Of course, the summation convention applies as usual to the Repeated 

( s) 

indices in the foregoing expression. Note that is an even function with 

(T) 

respect to the independent variable Ax. Likewise R^^ is even with respect to 
At. These important symmetry properties are a consequence of defining 
R^j according to Eq. (-6), rather than according to the alternative (-7) 
or (-8) . 

Because of their fundamental character and intimate relation to 
the turbulent energy, the scalar correlation functions and 

provide the natural bases for establishing appropriate length and time 
scales of the turbulence in the neighborhood of the space/ time point x,t. 

Specifically, let dv* = dxj^dx^dx^ represent an infinitesimal 
element of volume, and dt* an infinitesimal element of time at the variable 
point P*. Also let 

Ax* Ax * Ax^Ax^ = (AXj^)^ + (Ax^) ^ + (Ax^)^ 
represent the square of the spatial separation PP * . 

Then the desired characteristic length X* and time t* may be defined 
as follows : 



- A*^(x,t) 



/ 1 I (Ax*Ax)dv' 
|dv' 



T*^ = T*^(x,t) 



/lR^^^l(At)^dt' 

/iR^pIdt' 



(- 11 ) 



(- 12 ) 
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The integrals in (-11) extend over all space at a fixed time t* 

The integrals in (-12) extend over all time at a fixed spatial point x. 

Of course, these integrals are finite, despite the infinite domains of 
integration, because of the rapid decay in the correlation functions at 
large separations. 

Note that the integrands contain only the absolute values of the 

correlation functions; this ensures that every non-zero element of the 

correlation, even if locally negative in some places, makes a positive 

2 2 

contribution to the integrals which define X* and t* • 

The numerators of Eqs. (-11) and (-12) represent the second moments 

of the respective correlation functions. They are therefore analogous 

2 2 

to moments of inertia. Also, the quantities A* and t* are analogous to 

variances in probability theory. It is well known that such second 

moments attain their minimum values when evaluated with respect to the 

centroidal point of the distribution in question. In the present instance 

the correlation functions are symmetrical, with their centroids at the 

respective origins of coordinates, namely at Ax = 0 and at At = 0. Hence 

2 2 

the characteristic quantities A’^ and t* represent the minimum possible 
’’variances" of the respective correlation functions. These properties 
make A^^ and t* particularly appropriate measures of the length and 
time scales of turbulence in the vicinity of the arbitrary space/time 
point X, t . 



2.4 SPECTRUM OF TURBULENCE 



The space/time correlation tensor may be transformed into the 
wave number/ frequency domain by means of multi-dimensional Fourier 
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integrals. Let k\ <2 K 3 be the components of the spatial wave number vec- 
tor ic. Let 0 ) be temporal wave number, that is, angular frequency. Also 
let dv^ = dKj^dic^dK^ represent an infinitesimal element of ’’volume" in 
wave number space. Note therefore that 

lc*Ax = k^Ax^ = k^Ax^ + (2.4-1) 

Me may now define the Fourier transform pair 



$ (x,t;ic,iij) = 



„ f- - V +1(k*Ax + (jjAt)j I 
R^j(x,t;Ax,At)e dv dt 



(- 2 ) 



R^^(x,t;Ax,At) = [ 

/9'rT J J 



f 



. ^ - \ *“i(ic*Ax + o)At)j , 

. . (x,t;K,u))e dv dw 



jT, 

The integrations in Eq. (-2) extend over all of space/ time Ax, At 
with parameters x,t held constant. The integrations in (-3) extend over 
all of wave space with parameters x,t again held constant. 

The spectral correlation function like its Fourier transform 
is a nine component tensor. Naturally these two related functions 
exhibit corresponding symmetry properties. 

The first invarient 4> , . like its Fourier transform R.. is a scalar 

ii li 

and has a corresponding basic significance. It defines the character of 
the energy spectrum in the vicinity of the space time point x,t. 
Incidentally, it can be sViOwn from energy considerations that must be 
non-negative everywhere. 

It may be seen that the functions R^^ and are exactly equivalent 
in the information they provide concerning the structure of the turbulent 
field. The reason is simply that if either one of these functions be 
specified, the other one can be found from the above Fourier integral 
relations. 



(-3) 
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It is clear therefore that fundamental scales of length and time, 

(s^ 

which were previously developed in terms of and can equally 

well be defined in terms of the corresponding transforms and . 

Because of the above equivalence between the correlation functions 
R. . and the spectral functions it is immaterial in terms of which of 
these our turbulence theory is developed. We have chosen to emphasize 
the correlation aspects, and it is seen that there is no loss of infor- 
mation involved in so doing. 

It should be noted that the nine spectral functions (x, t ;Ax, At) 
are ostensibly eight dimensional, the independent variables being x^, 
x^, t, Ax^, Ax 2> Similarly, the nine spectral functions 

(x,t;<,o)) are also eight dimensional, the independent variables being 
x^, x^, x^, t, K^, K^, 

In addition to the above eight Independent dimensions which are 
explicit, there is an additional dimension which is implicit. Recall 
that the correlation functions are themselves averages theoretically 
computable from a large ensemble of individual flow realizations. We may 
assign to each realization a superscript index n which varies from one to 
some very large number. Then n is equivalent to an additional independent 
dimension. For if we wish to know not just the ensemble averages them- 
selves, but also something about their statistical distributions, then it 
becomes necessary to take into account the possible variations of the 
independent realization n from the ensemble mean. 

On the other hand, any specialization of the conditions reduces 
the number of independent variables which must be considered. Thus, for 
stationary turbulence, becomes independent of t (but not of At) For 
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completely homogeneous turbulence, becomes independent of x^, x^, x^. 

Also, for isotropic conditions, spherical symmetry prevails, and 
ceases to depend on Ax^ , Ax 2 » Ax^ individually, but becomes a function 
only of their resultant >/(AXj^) ^ + (Ax 2 )^ + (Ax^)^. Furthermore, in this 
case the number of distinct components of reduces from nine to two, 
depending only on whether the indices i,j are equal or unequal. These 
are vast simplifications. They account for the fact that homogeneous 
and isotropic turbulence has been studied so much more extensively 
than any other case. 

It should be noted that the fundamental phenomena of turbulence as 
represented by the ensemble of functions u^^\x,t)have five independent 
variables, namely, x^, ^ three dependent variables, namely, 

Uj^, u^, and u^- It is a curious fact that when we seek to make these 
phenomena more comprehensible by means of the correlation functions 
the number of dependent and independent variables, instead of being 
sharply reduced, are actually increased! Of course, the functions R^^ 
have a regular and symmetrical structure, they vary smoothly, and die 
away rapidly at large separations, whereas the u^ vary in an extremely 
complex and erratic way. It is this qualitative difference that makes the 
correlation functions seem so much more comprehensible. 



2.5 RELATION BETWEEN TURBLTENT CORRELATION SCALES AND MEAN FLOW 



It is customary and in many ways convenient to resolve fluid 
turbulence into two distinct complementary aspects which we then refer to 
as ”.the mean flow" and "the turbulent fluctuations", respectively. While 
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such a distinction is useful for certain conceptual and computational 
purposes, it should not be forgotten that these two aspects are inextri- 
cably coupled within a single unified process. In this connection it is 
pertinent to note that while we often find it convenient to introduce 
this distinction between mean flow and turbulent fluctuation into the basic 
Navier Stokes equations, the distinction is not innate within the original 
equations themselves. This is shown by the fact that these equations, 
being non-linear, cannot be fully decoupled into two independent sets, 
one dealing exclusively with the mean flow, the other dealing exclusively 
with the turbulent fluctuations. 

The foregoing fact provides a useful hint in connection with the 

problem of formulating appropriate auxiliary hypotheses, as required for 

supplementing the averaged equations of motion. We have seen that any 

such hypotheses will naturally be formulated in terms of the fundamental 

physical parameters which characterize the overall turbulence in the 

vicinity of a given space/time point x,t. Among these quantities are 

2 

the mean square velocity q (x,t) and the macroscales of turbulence in 
space and time, as expressed by the previously defined parameters X* and 
T* . 

Of course, in most problems of practical importance, interest centers 
chiefly on the mean flow field. One reason for this is that the mean 
flow is usually much more evident physically, much more convenient to 
observe and measure, seemingly more easy to understand. We tend to 
regard the turbulent fluctuations as a nuisance which has to be dealt with 
in order to get to the mean flow, and we tend to confine attention to 
those aspects of the turbulence whose connection to the mean flow is 
obvious and unavoidable. 
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More often than not, we do not have available the data which would 



be needed to establish such largely unseen characteristics of the turbul- 
ence as the spatial macroscale X* and temporal macroscale t*. Therefore, 
it is of doubtful utility to formulate any heuristic theory of turbulence 
directly in terms of these quantities. For example, what good does it 
do to express the Reynolds stresses, say, in terms of X* and t* if the 
latter quantities are initially just as much unknown as the Reynold's 
stresses themselves? Consequently, we must at this stage of our knowledge 
replace the fundamental macroscales X* and t* in our formulation by other 
quantities which are more or less equivalent, but which are far more 
accessible to observation and calculation. 

Fortunately, for this purpose we can take advantage of the intimate 
connection as noted above, which necessarily exists between the local 
turbulence field and the local mean flow field. Presumably because of 
this connection, these local scales of space and time must somehow be 
reflected in the local characteristics of the mean flow, and must 
therefore be deducible from the mean flow. 

While the exact nature of this relationship between the macroscales 
and the mean flow is far from clear we can legitimately infer that the 
relationship is inherently somewhat local in character. This follows from 
a fact universally revealed by all observed correlation curves. These 
always show a rapid decrease in correlation with increasing separation in 
space and time. Moreover, the correlations become negligible at separa- 
tions larger than about two or three multiples of X* or t* . 

Consequently, we infer, or at least we hypothesize, that quantities 
more or less equivalent to the true macro seales X* and t* can be defined 
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and computed in terms of the observable characteristics of the local mean 
flow. Let us call these quantities the apparent macroscales and designate 
them by the symbols A and i to distinguish them from the true macroscales 
A* and i* as earlier defined. The important problem of just hovr A and t 
should be defined in terms of the observable mean flow characteristics 
is taken up in a subsequent section. For the present we observe simply 
that if A and i be defined in some suitable manner, then it seems probable 
from the above reasoning that the dimensionless ratios 




will be very stable and well behaved universal functions which vary but 
slowly if at all. We assume that for many practical applications it will 
turn out to be premissible to approximate these universal functions by 
suitable average values, and to treat them as simple qua«i-constants . The 
actual numerical values of these ratios is regarded as a matter to be 
determined empirically. Of course, these numerical values will depend on 
the exact definitions that are finally adopted for A and t. However, the 
presumption is that if the definitions adopted are near optimal, the 
ratios will turn out to be roughly of order unity. 

It follows from the foregoing discussion that our definitions of the 
apparent macroscales A(x,t), i(x,t) in the vicinity of an arbitrary space/ 
time point x,t should depend only on mean flow conditions in a finite 
region surrounding that point. This principle is more realistic than any 
method which seeks to define these macroscales in terms of mean flow 
quantities only at the point x,t itself. On the other hand, it avoids 
the opposite extreme which would hold that A(x,t) and i(x,t) are somehow 
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directly dependent on conditions at all points in the flow field, however 
remote in space or time. Our proposed model accords more nearly than 
either of these alternatives with the observed fact that turbulent 
correlations are effectively limited to certain finite separations. 

Naturally, the definitions of X(x,t) and i(x,t) should emphasize 
mean flow conditions in the immediate vicinity of the point x,t, should 
give progressively less weight to conditions farther away, and should 
finally neglect mean flow effects at points sufficiently remote in space 
or time from point x,t. This requirement suggests the general idea of 
defining X and t in terms of appropriate weighted averages of certain 
pertinent mean flow quantities. 

The preferred choice of v/eighting functions for this purpose would 
be the correlation functions and themselves, if these were 

known. Of course. If these functions were knoi>7n, they could be used to 
find the true macroscales A* and i* directly, and there would then be no 
need for weighting functions. However, they are not known. Consequently, 
w^e are obliged to used heuristic substitutes, preferably functions which 
resemble the correlation functions to a certain extent. Fortunately, 
great accuracy is not necessary in this regard, for weighted averages 
tend to be relatively insensitive to minor variations in the form of the 
weighting functions used. 

In view of these considerations, vjc choose the Gaussian curve as 
an appropriate form of heuristic weighting function. It expresses the 
general trends of correlation in a suitable manner, and has various 
convenient and known mathematical properties as well. Consequently, we 
finally write the space/time weighting function as follows 
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-3 



Ax»Ax _ (At) 

.2 2 

W = e ^ ^ (-3) 

The reason for inserting the nuinerlcal factor 3 in the first term of the 
exponent will be explained presently. 

Of course, for purely spatial weighting, we may delete the temporal 
term in the exponent and set 

Ax* Ax 



-3 



= e 



Similarly, for purely temporal weighting, v;e may delete the spatial term 
in the exponent and set 

_ 

wW = e 

The above weighting functions and may be regarded as 



(-5) 



( s) 

initial heuristic approximations to the true correlation functions 
(T) 

and y respectively. Note that all four of these functions equal 

unity at zero separation. 

If we now replace the true correlation functions by their heuristic 
Gaussian approximations in Eqs. (2.3-11) and (2.3-12), we obtain the 
interesting results 

= A^ (-6) 

T*^ = (-7) 

This shows that for Gaussian correlation functions with exponential 
constants A and t, the true macroscales A* and t* may be simply identified 
with these exponential constants. Of course, this simple identification 
is possible only if the numerical factor 3 is inserted into the spatial 
exponent as indicated. This arises from the fact that the spatial integra- 
tion of Eq. (2.3-11) is three dimensional, whereas the temporal integration 
of Eq. (2.3-12) is only one dimensional. 
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Recall, however, that the actual correlation functions are not 
necessarily Gaussian and are in fact unknov:ns in the present analysis. 
Moreover, for the reasons previously given, the quantities X and t will 
be defined not in terms of the correlation functions, but heuristically 
in terms of the mean flow. In fact the basic reason for introducing the 
above weighting functions in the first place is that they are needed in 
connection with these heuristic definitions of X and i. However, the 
weighting functions in turn contain X and i as parameters. It is there- 
fore apparent that the actual calculation of A and t v^ill necessarily 
involve some type of iteration process. 

For certain purposes in the subsequent analysis it will prove 
convenient to utilize a normalized form of the above Gaussian w^eighting 
function. In this connection, consider the following integrals, namely. 



H = /AMv*dt’ (-8) 

= //W^®^dv'dt' (-9) 

and 

= //W^'^^dv’dt (-10) 

As usual, the Integrands extend over the entire space/ time domain. Of 

course, the integrals themselves are alv/ays finite even if the domain of 
integration extends to infinity. 

We now define the normalized weighting functions as follows: 



VJ 



W 

H 



W 

V — 

H 

(T) W 
w = — 

II 



(s) 

(T) 



(- 11 ) 

(- 12 ) 

(-13) 
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When normalized in this way, the above weighting functions have the 
following convenient properties, that is, 

//wdv’dt’ = 1 (-1A) 



//w^®^dv'dt' = 1 
//w^’^^dv'dt' = 1 



(-15) 

(-16) 



2.6 LOCAL LENGTH AND TIME SCALES OF THE MEAN FLOW 



We have been considering any arbitrary and general turbulent flow 
field of which the mean flow may be either steady or unsteady. With every 
space/time point of this field may be associated two parameters X and i 
which express respectively a length scale and a time scale characteristic 
of the mean flow pattern in the vicinity of the point. Our problem is to 
devise suitable explicit definitions of these parameters in a way which 
meets various essential requirements. Some of these have already been 
discussed. There are some additional requirements as well, among which 
are the following: 

(1) The parameters X and t must be true scalars, and therefore invarient 
with respect to any rotation, reflection, translation or acceleration 
of the reference axes. 

(2) These parameters should also be everywhere continuous, finite and 
positive (except possibly in certain limited singular regions, such 
as at a solid boundary) . 

2 2 

(3) If feasible, the mean flow quantities X and i should preferably be 

2 2 

related to the correlation quantities X* and t* , respectively, such 
2 2 

that the ratios r and — ^ are close to unity, or at least so that 

X* T* 

these ratios remain as nearly constant as possible. 
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These requirements can be satisfied as follows. Let 



Uj = U^(Xj^,X2,x^,t) = u^(x,t) 



( 2 . 6 - 1 ) 



represent the velocity field of the mean flow. Then the corresponding 
strain rate tensor is 

fdV. 



j ■ 



J 

Of course, we are considering here only the incompressible case, 
for which the first invariant vanishes, that is 



Til = 21 





- 2rf^i 


1 




1 




1 




Li®”!. 




3 x» 

2 J 




[3x3] 





0 



(- 2 ) 



(-3) 



Consequently, there is no net dilatation, and the represent purely 
distortional effects. 

We now define a generalized strain rate Qy and a generalized strain 
rate gradient by means of the following expressions, namely 






.2 




(-5) 



From these definitions, the following relation can be obtained. It 
will prove useful in the subsequent work. Thus 




The factor 1/2 which appears in Eq . (-4) calls for a brief explanation. 
Consider a simple plane parallel shear flow such as illustrated in Fig. 

2.6-1. Let the velocity be parallel to axis x^, and let the shear gradient 

be in the direction of x^* Then F^^^ = T = F^^ “ ^^3 * ^ 3 ^ " 

the only non-zero shear strain rate is F^^ " ^^l' applying Eq . (-4) 
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y 




In this case the definitions of and 
Q,' reduce to 

XI = (dU/dy) 

Xi'= U" = (d^U/dy^) 



FIGURE 2.6-1 EXAMPLE: SHEAR QUANTITIES XI AND X2' 

IN PLANE PARALLEL FLOW 
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we find that, because of the presence of the factor 1/2, we obtain the 
simple result 



In other words the presence of the optional factor 1/2 identifies Q 
with the ordinary shearing rate in a simple plane shear flov. Eq . (-4) 
may be regarded as an extension of this comcept to a general state of 
strain rate. 

It can be shown that while U, and f,. depend or the orientation of 

i ij 

the reference axes, the quantity Q is invariant in this regard, and is 

12 2 

therefore a true scalar. The same is true of the quantities 0. y Q ^ Q' 

2 2 2 

and . Hence X and t can be convenient!}^ defined in terms of any 

independent pair of these quantities. 

To illustrate this idea, consider the expressions 



X 



2 




(- 8 ) 



2 

T 



>2 



(-9) 



2 2 

v/here and C:^ are dimensionless constants. The combinations shown may 

be seen to define local scales of length and time, X and t. In fact, 

Eq. (-8) v^ill be recognized as the expression for the von Karman mixing 

length, or rather, as a generalization of the concept, for von Karman's 

original definition v^as restricted to the special case of parallel flow. 

One difficulty with the above approach, however, is that the functions 
2 2 

A and t so defined tend to fluctuate too wildly instead of being smooth 

2 2 

and well behaved. Moreover, at certain singular points where Q or Q' 

2 2 

are locally zero, the quantities X or t may reach values of zero or 
infinity, neither of which is premissible or physically meaningful. 
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The cure for this difficulty, of course, is to replace the local 
2 2 

point values of Q and by suitable space/time averages evaluated in 
the vicinity of the point. Thus we may write 



= //wfi^dv'dt' 



(- 10 ) 



0^^ = //wQ '^dv’dt' 



(- 11 ) 



whereupon 



X 



2 



,2 9. 






(- 12 ) 



2 

T 



= c 



7 



(-13) 



The nevr values of X and x obtained in this way will now vary 

smoothly as required, and will never be either zero or infinite. The 

2 2 

dimensionless constants and remain to be specified. 

Note that the normalized form of the weighting function is used in 
the integrals (~10) and (-11) . This is important because it ensures that 
the integrals remain relatively insensitive to the precise form of the 
weighting function . 

For the reasons given earlier, a Gaussian curve has been chosen for 
the weighting function. For consistency, the characteristic exponential 
constants X and x of this function should be respectively equal to the 
X and X computed from Eqs. (-12) and (-13). In other words, the 
averaging intervals in length and time should be approximately equal to 
the characteristic local scales of length and time. Some iteration will 
be involved in satisfying this requirement. 

It turns out to be useful to generalize the above idea even further. 
For this purpose we define the following integrals, namely. 
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(-1M 



= //w(fi^)‘^(fi'^)‘^dv’dt' (-15) 

where the exponents a, b, c, d may be assigned arbitrarily. However, 
for reasons which will soon be evident, we stipulate that 



D = 



b 

(a -h b) 



d 

(c -h d) 



¥ 0 



(-16) 



Moreover, to avoid improper integrals we require all of the above exponents 
to be non-negative. 

Now, depending only on how the arbitrary exponents a, b, c, d are 

2 2 

specified, the integrals I and J may be assigned various physical 

dimensions. Consequently, it is always possible to define a mean square 
2 2 

length A and mean square time x by expressions of the following general 
form, namely. 



CM 


= C^(I^)''(J^)^ 


(-17) 


2 

T 


= C^(I^)"'(J^)" 


(-18) 


2 

where and 


2 

are arbitrary dimensionless constants. 


It is easily shown 



on dimensional grounds that the four exponents k, i, m, n are given by 
the following formulas: 



k = --i(c + d) m = +^ 

(-19) 

i = + ^(a + b) " = - D 



It is now seen that Eq . (-16) represents a necessary condition to 
ensure the existence of the above solution. 
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The foregoing formulas show that, subject only to the condition (-16), 

2 2 

specification of the six arbitrary constants a, b, c, d, C. and C uniquely 

A T 

2 2 

defines the quantities A and t • These six ouantities are here regarded 

as empirical constants whose values can in principle be chosen so as to 

2 2 

confer upon the defined quantities A and t a wide range of desired 
properties . 

Ideally, it would be desirable to choose these six constants in 

2 2 2 
such a way as to make A and x agree as nearly as possible with A* and 

2 

X* , respectively, over a broad range of conditions. In principle, this 
can always be done if and when sufficient experimental data becomes 
available for this purpose. At the present time, the available data are 
insufficient for such an analysis. It therefore becomes necessary, at 
least for the time being, to resort to expedients of a more heuristic 
nature, as will be explained presently. 

2 2 2 

The calculation of the space/time integrals I and J which fix A 

2 

and X is seen to be inherently lengthy and complicated. One reason is 

that the space/ time integrals involved are, in their present form, four 

dimensional. The second reason is that these same quantities appear not 

only as the final unknowns, but also as the exponential constants in the 

Gaussian weighting function. Hence our next endeavor is to simplify the 

analysis as much as possible in both these respects. 

In this connection, we observe first that in computing the integrals 
2 2 

I and J , if the mean flow happens to be steady, the integration over 
time cancels identically from the result, and the integrals reduce to 
ordinary volume integrals. Secondly, we notice that if the flow field also 
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happens to be either plane or axi-symmetric, all volume integrals can be 
reduced to corresponding surface integrals. Also if the flow field happens 
to be of the type known as self-similar , all volume integrals again 
reduce to surface integrals. If the flow field is both axi-symmetric and 
self-similar, then a double reduction is possible, and the volume integrals 
can be reduced to line Integrals. These are tremendous simplifications. 

By a self -similar flow field, we mean a field in which there exists 
a family of surfaces across which the mean velocity distribution has a 
common pattern; the velocity distribution across the various surfaces 
differs, at most, only in the respective scales of length and time, or 
what amounts to the same thing, of length and velocity. Many common flows 
in jets, wakes, ducts and boundary layers are known to be self-similar. 

The characteristic surfaces in such cases are conveniently taken as cross- 
sections approximately normal to the principal flow direction. 

This consideration of self-similarity shows us, thirdly, that for 
such flows each characteristic surface possesses its own overall length 
and time scales; these are constant over, any one such surface, but may 
differ from one surface to the next. In addition to these constant 
overall scales of length and time which characterize the cross-section as 
a whole, there are also distinct local scales of length and time, as 
considered earlier, in the vicinity of any reference point on the surface. 
These local scales, unlike the overall scales, vary smoothly from point to 
point along the cross-sectional surface. 

Of course, even when we have taken these factors into account to 

2 2 

simplify our formulas, the definitions of X and t will remain indeter- 
minate until the four exponents a, b, c, d and the two proportionality 
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constants C. and C are definitely specified. Obviously, there exists an 
enormous amount of flexibility here, and this was done intentionally. 

The flexibility has been deliberately provided in order to permit matching 

2 2 2 2 

A and t as closely as possible to A* and t* over a broad range of 

conditions . 

Naturally, to do this properly will require much data beyond what is 
now available. Hence at the present time our scheme has far more 
flexibility than we really need or can properly exploit. 

Consequently, we have found it necessary to simplify at this time in 
what may seem to be a somewhat arbitrary manner. Actually, such simplifi- 
cations were guided by a fair amount of numerical experimentation. Accord- 
ing to present evidence, they seem to give reasonably satisfactory agree- 
ment with physical experiment. No doubt, there remains room for further 
improvement, however. 

Specifically, for our present turbulence model, we have chosen the 
following numerical values for the four arbitrary exponents, namely, 

a=2 b=0 c=l d=l (-20) 

Also, it proves to be most convenient to choose the proportionality 
^ factor simply as 

cj = 1 (-21) 

2 

The other proportionality factor is not required for the present simpli- 
fied analysis. 

Let us summarize the situation as it now stands after the foregoing 
simplifications. The integrals below extend over the entire flow field. 

Thus 

Ax = (x’ - x) = separation variable (“25) 



71 






W(x,3c*) = e = V7eighting function 

H(x) = /Wdv* = normalizing integral 
w(i,x) = -^ W = normalized weighting function 
2 4 

I (x) = /wf2 dv* = first characteristic integral 
2 ^ 2 

J (x) = /w(f2Q*) dv* = second characteristic integral 

9 

2 1 “ 

X (x) = — 2 = local length scale parameter 



(-26) 

(-27) 

(-28) 



(-29) 



(-30) 

(-31) 



Isotice that Eq. (-26) above differs from Ea . (2.5-4) in that the 
unessential factor 3 has been eliminated from the exponent for convenience. 
This is arbitrary but permissible; it amounts to a rescaling of X. 

The foregoing relations suffice in principle to define the required 
scale factor X at every point in the flow field as required, but do involve 
a computational difficulty. The problem is that the value of X Is needed 
for Eq. (-26) before a computed value is available from Ea . (-31). Hence, 
it is necessary to resort to an iteration process such as the following. 

Let the (n 1) th approximation to X be defined as: 






(-32) 



where 



(li) = / 11 W 

n+1 ‘'all space n 



(x,3?')J2^(3^')dv' 



J . 1 (x) = / W (x,x*) [Qfc * (x’) ] dv’ 

n*fl •'all space n ’ 



The weighting function is defined as follows. Let 

W (5^.5^.) = exp 



. (-33) 
(-34) 



(-35) 
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Then 




n 



w (x,x) = 



n 



W (x,x*)dv* 



(-36) 



/ 11 ^ 

*' all space n 



We nox'7 define (in principle) 




(-37) 



n-x» 



Advantage should be taken, of course, of any Inowledge concerning the 
flov to hasten the convergence of the X *s. In particular, if the flow is 
self-similar, such as the flow in a pipe, duct, or boundary layer, or the 
far downstream region of a turbulent jet or wake, calculations need only 
he perform.ed on a single representative cross-section of the flow, and 
the results rescaled according to the appropriate nondimensionalizing 
parameters . 

As an initial guess to start the iteration process, we take 
everyw’here. Then in computing X^, we find that the weighting function 
is equal to unity everyv^here, and so we obtain simply 



approximation to X that we obtain. 

Numerical experiments with various flow fields have shown that, in 



in Table 2.6-1 vjhich shows successive numerical approximations to the X 
distribution in laminar flow in a two dimensional channel. Consequently, 
for the purposes of the calculations presented in this paper, the 
approximation was made: 




4 



(-38) 



a constant independent of position. Thus, ^2 first nonconstant 



general, the convergence of the is extremely rapid. This is illustrated 



X(x) X^Cx) 



(-39) 
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Some improvement could possibly be gained by continuing to higher 
approximations, but in view of the gross simplifications already inherent 
in the model as a whole, the authors feel that such refinement is not 
warranted, ^1oreover, absolute accuracy in this regard is not essential. 
It is necessary, however, that the chosen order of approximation be 
adhered to consistently. It is also important to adjust the initially 
undetermined numerical constants of the model so that the chosen order of 
approximation gives satisfactory agreement with experimental data. The 
significant numerical constants are those that occur in the heuristic 
expressions for the quanities a, 3, and y- 



TABLE 2.6-1 

Convergence of Macroscale Distribution in a Two-Dimensional Channel 
with Parabolic (Laminar) Velocity Profile 



y/b 






A /b 
n 








n = 0 


n = 1 


n = 2 


II 


n = 20 


0.0 


oo 


0.775 


0.695 


0.674 


0.664 


0.2 


00 


0.775 


0.706 


0.693 


0.689 


0.4 


00 


0.775 


0.731 


0.727 


0.727 


0.6 


00 


0.775 


0.758 


0.758 


0.758 


0.8 


00 


0.775 


0.782 


0.781 


0.781 


1.0 


00 


0.775 


0.802 


0.800 


0.800 


^ 2b = 


channel v;idth 


y = distance 


from centerline. 
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2.7 REYNOLDS STRESSES MD GENERALIZED EDDY VISCOSITY CONCEPT 



In dealinp with an incompressible fluid, it is customary to divide 
the equations of motion through by the fluid density. As a result, the 
viscosity, pressure and stress terms all appear in these equations divided 
by the density. It is convenient to avoid siting out the density explicitly 
each time. Therefore we define the quotient of viscosity divided by 
density as the kinematic viscosity, and assign it the symbol v. Similarly, 
the quotient of pressure divided by density becomes kinematic pressure, 
symbol p. Stress divided Vy density becomes kinematic stress, symbol 
Note that kinematic pressure and kinematic stress have dimensions of 
velocity squared. 

The kinematic Reynolds stress tensor may now be written 



T . . = T , . = “U ,u . 
13 i 3 



(2.7-1) 



The contraction obtained by equating indices and summing over the 
repeated index defines a scalar invarient, namely, 

'ij ■ 'll + '22 + '33 ■ -'“i + “2 + “3’ ■ ■ ■ 2 '= 

How the Reynolds stress tensor may be rewritten as the sum of a 
purely dilatational tensor plus a purely distortional tensor, that is, 



(- 2 ) 



T.. = - -:pq 6 .. + T.. 

13 3 13 13 

Conversely, the distortional part is defined by 

T . . = T . . + -rq 6 

13 13 3 ‘ ij 



(-3) 



(-4) 



Of course, the individual components of depend on the orientation 



of the coordinate exes. However, the double summation 

1.1 »2 

2 13 13 



(-5) 
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is invarient with respect to any change in the orientation of the reference 
axes at any point, and is therefore a true scalar. We may think of i’ 
as representing a kind of generalized scalar amplitude of the tensor 

Of course, the optional factor 1/2 was inserted into Eq . (-5) to 
maintain consistency with the corresponding formula for the strain rate 
invarient The physical justification is as follows. Consider a simple 
state of pure shear stress such that = 1^2 = '^^93 = 

simple shear 1^2 ” '^21 only non^zero stress. In this case Eq . 

(^5) gives 




Hence the invarient i’ may be identified with the magnitude of the simple 
shear stress itself. Eq . (-5) represents an extension of this concept to 
the general state of stress. 

A somewhat similar analysis may be made of the rate of strain tensor 
of the mean flow field. Recall the definition 



r 



ij 




(-7) 



For an incompressible fluid, the associated dilational tensor vanishes, 



that is. 



^il ^11 ^22 ^33 



rSU 8U 3U« 

^ + T— ^ ^ 



9x, 



9x. 



9x, 



= 0 



(- 8 ) 



Consequently the tensor is itself wholly distortional , 



An amplitude 



scalar U may be defined through the invariant relation 



1 

2 



r r 

ij Ij 




(-9) 
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Next we wish to formulate an appropriate way to express the general 



relation that may exist between the distortional stress tensor and the 

distortional strain rate tensor F,.. In this connection, it is instructive 

ij 

to consider as a preliminary model, the relation which is known to exist 
between the ensemble averaged distortional viscous stress 
distortional strain rates F^^ of the mean flow. This is given by 



(t! .) = VF. . 

ij V ij 

where v is the ordinary molecular kinematic viscosity. 

This model suggests that we attem.pt to express the distortional 
Reynolds stresses in the form 



(- 10 ) 



t: . = ef. . (-11) 

IJ IJ 

where v is now replaced by e, the so-called eddy kinematic viscosity. 

However, this form of equation implies that 



T * T * T * 

11 _ In ^ 33 

F ” F F 

11 22 33 



T ^ T * T ^ 

^23 ^ ^ ^ 12 

TTY 
23 31 12 



(- 12 ) 



Unfortunately, there is no a-priori reason to believe that the tensors 
and F^^ will necessarily satisfy the five similarity constraints 
expressed by Eqs. (-12). We therefore require a more general relationship, 



one which remains applicable even if and F.. are both specified 



iJ 



IJ 



arbitratily . 

One method which is sometimes suggested for dealing with this 
requirement is to postulate a relation of the form 



T ’ = r r 

^Ij ^ijkl kl 

Now the eddy viscosity becomes a fourth order tensor. Since t 1 , and F.. 

^ ^ ij ij 

each have six distince components, the tensor e , has thirty-six 



(- 13 ) 
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components! By employing tensor arguments, the number of independent 
eddy viscosity components can actually be reduced to ten. Furthermore, if 
the tensors and are arbitrarily specified, Eqs. (^13) can in fact 
be reduced to an equivalent set of five independent equations in ten 
unknowns. However, these equations do not constitute a determinate set. 
Clearly, this type of formulation is both inadequate and overly complex 
as well . 



A somevhat simpler tensor relation has been suggested by Hirt, 
Ref. (A). In our notation it becomes 



’Ij ' 'jk''kl’ 



(-14) 



While this is an Improvement over Eq. (-13), it is not free of difficulties. 

If is taken as an unsymmetrlcal tensor, Fq . (-14) falls to produce a 

detenninate solution for e, ^ when xl. and F , . are specified. On the other 

kl ij Ij 

hand, if is taken as symmetrical, the relation (-14) loses generality; 

there are then certain theoretical combinations of x’ and F.. which are 

ij ij 

not reconcilable with (-14) no matter how the e , are chosen. 

J ^ 

The foregoing difficulties can be overcome by adopting the following 

relation between x* and F.., namely, 
ij ij 






(- 15 ) 



It will be shown that the quantities e and f^^ suffice to define the 
relation between distortional stress and strain rate in a fully general 
way. If x|j and F^^ are specified arbitrarily, a corresponding unique 
solution can be found for e and • Conversely, if e, 

specified, then j uniquely determined. 

Solving Eq. (-15) for f^^ gives 






(- 16 ) 
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It is immediately apparent from this that 
f., = f. 



ii 



ji 



and 



(- 17 ) 



(-18) 



f . . = f ^ T + f + f = 0 
11 11 22 33 

However, as it stands, Eq . (-16) does not uniquely fix f^^ because 
it contains the as yet undetermined parameter e. It will now be shown 
that for any arbitrarily prescribed and a corresponding "preferred" 

value of e can be uniquely defined. For this purpose let 



t'' . = cfif . . = T 1 . - ef . . 
ij iJ ij iJ 



Here represents the discrepancy between and 



= t’ 



t! . - 2 er, .t! . + e r, .f 



ij ij 



ij 1.1 



ij ij 



or 



7 ? 7 2 

t"“ = t' - er . .t' + e 



ij ij 



„2 



Now to minimize the invarient t of the discrepancy, we set 

...2 



3t 



dE 



= 0 - r . .T^ . + 2ef2 =0 
ij ij 



whereupon 



r ,. T ’ r ,. T '.. r ,. T .. 

= ij iJ = ij 11 = _ij_JJ. 

2fi ij ij ij ij 



(-19) 



We wish to 

choose e so as to Tninimize this discrepancy in some appropriate sense. 
Hence we form the invarient of the discrepancy namely, 

t’.'.t" = 2t"^ = [t! . - er, .Hi' . - eP, .] 

ij ij ij ij iJ ij 



(- 20 ) 



(- 21 ) 



(- 22 ) 



This simple result is important. It is the basic definition of 
the scalar eddy viscosity e for the case of general stress and strain rate. 

Now multiplying (^19) be and using (-22) we obtain the significant 



consequence that 



This is the mathematical expression of the concept that the tensors T 

ij 

and are ‘orthogonal*'. Hence our basic s tress/strain rate law, Eq . 

(-15), expresses the idea that the Reynolds stress tensor consists of 
a tensor ^hich is simply proportional to ^ plus a tensor 

which is orthogonal to 

The tensor f^^ has nine possible components. The symmetry constraint 

(-17) reduces this to six distinct components. The non-divergence constraint 

(-18) reduces it further to five independent components. Finally, the 

orthogonality constraint (-23) reduces the degrees of freedom of f^^ to just 

four. Hence e and f^^ together possess five degrees of freedom. This is 

precisely the necessary and sufficient number to express the general 

relation between distortional stress and mean flow strain rate. If we add 
2 

as a sixth degree of freedom, this then allows us to define the general 
relation btween overall Reynolds stress and mean flow strain rate. 

Since the foregoing equations are pure tensor relations, they remain 
valid for all possible orientations of the coordinate axes. It is particul- 
arly convenient to choose the principal axes of strain rate in order to 
bring out the essential concepts in the simplest possible form. We use 
asterisks to designate these axes. Also, there i s no loss in generality 
in ordering the subscripts so that ^ ^ ^ 33 - course “ ^21 ~ 

r *2 " Moreover, it happens to be useful to introduce the following 
auxiliary terminology. Let 



rjj - ^ sin |I + 9* 



^ 22 " 71 ^ S±n(0 + e*) 



(-2A) 



r*3 - ^0 sin(-|l+ e. 
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The merit of Eqs. (-24) lies in the fact that they express the three 



quantities ^22’ ^33 terms of just tv^o parameters, namely, U and 0* 

Furthermore, as the reader may easily verify, eqs. (-24) satisfy identically 

the following important constraints, 

rye -f + r* =0 
^11 ^22 ^33 



1/2 2 2 

±\ pyc^ + pyc^ + pyc^ 

2 ' 11 22 33 



= Q 



(-25) 



The parameter G* is then fixed from the relation 



7 3 r*2 



= tan e* 



(r* 

^ 11 



r*3) 



(-26) 



It is equally useful to adopt a similar scheme for the • Thus 






f* 


2 


22 


7 3 


f* 

33 


2 
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g* Sin 



-|1+ 



(-27) 



wbareupon 



f* + f* + f* =0 

11 22 33 



(-28) 



o o o 

f*^ + 

11 22 33 



= g* 



and 



^22 
‘'h - 



= tan tp* 



(-29) 
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Now the orthogonality constraint, en . (~23) , reduces In this case to 



the form 



r* f* 4- r* f* -h r* f* =0 

11 11 22 22 33 33 



Upon substituting Eos. (-24) and (-27) into (-30) and reducing, we 
obtain the beautifully simple result 

If)* = 0* ± Y 

This too expresses a kind of orthogonality between the T and the 

With ip* now known, the components > ^22’ ^33 fixed as 

soon as g* is specified. The off-diagonal terms f*^* ^12 

specified individually. The final stress/s train rate equations for these 
particular axes now taVen the form 

= “fill 

T*„ = -u^ = - 4 + e[r*2 + 



22 



3 

1 2 



22' 



T-* = -n"" = - Cl + \ri( + Of* 1 

33 3 3 ” ^ 33 33^ 



^23 " ~^ 2“3 “ ® ^^ 23 ^ 



^31 ~ ~“3^1 “ ® 



t*2 = ® ^^31^ 



Note that the principal axes of stress and strain w^lll coincide if, 
and only if, f*^ = f*^ = 

The further reduction of these equations for the important special 
case of a plane mean flow is interesting. Now the axes x^ and x^ are in 
plane of the flow, and x^ is normal to this plane. Hence we have 
0 * = 0 



r*i = +n 



(- 30 ) 



(-31) 



(-32) 



the 



(-33) 
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0 



s= 

22 



= -r 
^ 33 



Consecuen tly from K*q . (-31) we find that 



= ± ~ 



(- 33 ) 



(-34) 



There is some choice concernin^^ the manner of handling alpehraic 
signs. Vj*e choose to drop the negative sign in (-34) and instead, to allow 
in the following expressions to tahe on positive or negative values as 



renuirecl- Then 






f* = - -7=i r* 

11 n 


f* = 0 
23 




fK =4- -5=5^ r'x 

"22 ^ 


f*^ = arbitrary 


(-35) 


f* = - -ri e* 
33 73 ^ 


0 

II 

CM 

-K«-H 





Of course, f*^ = f*^ = 0 by reason of symmetry. Hence g* and f*^ 
remain as the two arbitrary degrees of freedom which define the coefficients 
f*. . 

The final s tress/s train rate equations become 



T*i = -u^ = - ^ 0^ + eiUl 1 P,*] 

^22 “ ~^2 ^ ~ J 

"33 = -"*3 = " i ■ 4 



T 



* 

31 



-u^u^ = 0 + 



T 



* 

23 




0 



(-36) 



Hence the overall stress/s train rate relation is now fully specified 
1 2 

by the four quantities q , e, g*, and 
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In the case of axi-symmetric mean flows, the results are similar to 



those for plane flow, except that 0* is not in jxeneral zero; therefore ij>* 
must be found from Eq • (-31) and the f*^ from Eq . (-27). Hovrever, we still 



have f *2 = ^22 ” ^ this case. 



1 



The basic stress/strain parameters are now — q , e, g*, 

Note from Eqs. (-2Q and (-31) that ij;* is a known function of the strain 
tensor, but that the other four parameters may be arbitrary. 

For purposes of approximation in connection with our heuristic 
turbulence model, v/e have decided to make the simplifying assumption f^^ = 0. 
Actually this is not as drastic as it seems. In the first place, as we 
have seen, f^^ embodies at most only four degrees of freedom, not six or 
nine. Secondly, we are mainly concerned with plane or axisymmetric mean 
flows, for this f^^ has only two degrees of freedom, namely, g* and 
Thus this assumption amounts to neglecting but two degrees of freedom. 

Secondly, an analysis of the momentum transfer involved in generating 
the Reynolds stresses suggests that the term for most of the 

effect, and that the deviation . should indeed be small. While the 

ij 

experimental data w’hich is available is not decisive on this point, it does 
tend to support the above hypothesis. It does seem that this simplification 
may becom.e somewhat Inadequate in the immediate vicinity of a fixed wall, 
but it appears to be reasonably satisfactory elsewhere in the flow. 

Under these assumptions, the principal axes of stress and strain 
coincide, and the equations for these particular axes reduce to the simple 
form 



* = 
11 


1 

" 3 


2 

q 


+ cQ 


* s 

22 


1 

3 


2 

q 


+ 0 


* s 

33 


1 

3 


2 

q 


- C$7 


* s 

23 


X* 

31 




T* = 
12 



(- 37 ) 



84 



The foregoing stress/s train rate law, despite Its apparent simplicity, 
can be useful for calculating the Reynolds stresses which correspond to 
any specified mean flow strain rates only if the value of e can somehov^ be 
estimated with acceptable accuracy. Whether this can be done Is an 
empirical question, capable of being answered by suitable experimental 
information. The experimental data so far available Is not yet adeauate 
to answer this need fully. Further experiment and analysis is needed. 
Tentatively, however, we are justified in postulating that g Is some function 
of the characteristic local turbulence parameters q, X, and i as prevlsouly 
defined, perhaps also of molecular kinematic viscosity v, distance y from 
the nearest wall, and possibly of some as yet unknown additional parameters 
as well. This can be written symbolically in the form 

t = c(q, X, T, V, y, ) (-38) 

Actually, we would be on slightly more fundamental ground to choose 
X* and T* in place of X and t. Unfortunately, and t* are not usually 
knovm, v’hereas X and t have been related in a specified manner to the mean 
flow’. Conseouently , for pragmatic reasons we replace X* and i* by X and i 
In the above expression. 

The foregoing empirical law’ which governs the eddy viscosity e can be 
greatly simplified by dimensional analysis. We recall tliat in mechanics 
generally, the number of fundamental dlmesnions Js three. These are often 
taken as mass, length and time, llow’ever, in the case of Incompressible 
flow, dlvidion of the equations of motion by density has the effect of 
eliminating mass from the problem. Our final equations are therefore all 
kinematical in character, with length and time as the only fundamental 
dimensions. It follows therefore that any two of the above parameters may 
be chosen as the necessary and sufficient basis for the non-dlmenslonallzlng 
process . 
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From a purely logical standpoint, it would seem natural and convenient 
to use the basic local scales of length and time, A and x, for this purpose. 
However, there are grounds for believing that a and A play roles which are 
physically more important than that played by x. This seems to be especially 
true for mean flows which are steady, or which vary only slowly. Conse- 
quently, we choose q and A and our non-dimensionali zing parameters. 

As a result of the foregoing considerations, the empirical eddy 
viscosity law may now be rev:ritten in the dimensionless form. 






= a = a 



ill z 

X ’ V ’ X’’*' 



(-39) 



Thus, our problem has been simplified to that of finding the dimen- 
sionless coefficient a which presumably varies as a function of the dimen- 



sionless parameters 



ail and lai 



and 



, and possibly also of one or two 



other dimensionless variables as well. 

In practice, the usefulness of the above formulation hinges on whether 
the dimensionless coefficient a does in fact behave in a reasonably simple 
and stable manner. Our present indications are that this is indeed so. 

Owing to limitations, of time and data, we have not fully established the 
actual dependence, if any, of a on the quantities and . We have 

found heuristically , however, that satisfactory overall results can 
apparently be found by assuming that a depends only on the wall distance 
parameter (y/A) . The following expression seems to give reasonably good 
curve fits, namely. 



a = 0.065 



{. . .-|5 - ■!■) 



(-40) 



Of course, this expression is tentative. More data is needed on this 
point in any case. 
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Naturally, this particular expression for a is applicable only so 
long as we specifically retain our present definition of A, Any modifica- 
tion of the definition of A clearly calls for a corresponding adjustment 
in a. Note however that the dimensionless quantity 



( > 
e 


[ax' 


X*q 


[x* 



should remain substantially invariant with respect to such details. 

However, even if a be accurately known, and A as V7ell, the actual 

eddy viscosity e cannot be evaluated until q also is established. Bearing 

1 2 

in mind that the turbulent energy is E ci , 1 1 is apparent that the 
determination of q must involve an energy analysis of the turbulence. 
Physically, this is fundamental. It ties together the turbulence, the 
eddy viscosity and the mean flow into a consistent overall picture. It 
also provides the additional equations which are necessary to form a 
determinate set. 



2.8 HEURISTIC ENERGY /JS7.LYSIS OF THE TURBULENT FIELD 



Starting with the Navier Stokes equations of motion, it is possible 
to develop a corresponding exact energy equation, as shown in detail in 
Appendix A. The following result is obtained, namely. 




II 



III 



IV 



.r. . 

2 ij IJ 






9x . i 
1 



[U.E] - 



3x , 





u .u 

n + J J 




2 )J 



V 



( 2 . 8 - 1 ) 



+ 
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f^T. 


’'IJl 




3x^ 




3x . 
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This equation governs the evolution over time and distribution over 
space of the ensemble mean turbulent kinetic energy, that is, of 

r 1 2 1 1 , T . '2, , 

E = 2 'I =2 - 2 ^2 

The five terms on the right side of Eq . (-1) have the following 
physical significance : 

Term I - The rate of energy input into turbulence as work done by the mean 
flow against the resistance of the Reynolds stresses. 

Term II - The rate of irreversible energy dissipation into heat. Note that 
the quantity is positive definite. This ouantity accounts 

for the local entropy production. 

Term III-The rate of convective transport of turbulent energy E by the mean 

flow U . , 

1 

Term IV -The rate of net transport by turbulent diffusion of the turbulent 
pressure plus kinetic energy fluctuations. 

Term V - The rate of net transport of energy by molecular diffusion. Note 
that the Reynolds stresses, as well as the energy, appear in this 
term. 

While Eq , (-1) is theoretically exact, it contains several terms 
which cannot be evaluated exactly v/ithout detailed knowledge of the 
turbulent fluctuations. Such detailed Information is not practically 
attainable. It is therefore necessary to approximate these terms by suit- 
able semi-empirical expressions. These expressions vill necessarily contain 
empirical coefficients whose values can in principle be chosen so as to 
optimize agreement with available experimental data. 

We have shown in the preceding section that the Reynolds stresses 
can always be expressed with complete exactness and generality in the form 
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-u . u . 
1 J 



(- 3 ) 



where 



and 




1 

3 




+ e[r. . + Qf . .] 

iJ 



T. .r. 



c = 






r, .r. . 

ij iJ 



r . .f . . = 0 



(- 4 ) 



(- 5 ) 



Consecuently , term I, the l^crk input expression in the energy equation 
reduces exactly to the result 



1 



T . . r . . 

ij 



E cr^ 



(-0 



liov;ever, as \:e have seen earlier, e must here te approximated by an 
expression of the form 

e = a>q (-7) 

v/here a is presumably a well beliaved empirical function, the expression 
for which has been given earlier. 

Next consider the dissipation of energy into heat, as given by term II 
of tlie basic energy equation. Dimensional considerations suggest that 
this term can be expressed to good advantage in the form 



f:,. = 



vy • . Y . . 

IJ 




(-P) 



This expression amounts to a definition of the so-called dissipation length 
v:hich is a length characteristic of the microscale of turbulence. This 
form of expression would be very convenient if the dissipation length 
happened to be related to the macroscopic scale X in some fairly simple and 
invarient way. Unfortunately, this does not seem to be the case. 
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However, experimental studies such of those by Laufer, Pef. (^) , 
clearly suggest that at sufficiently high Reynolds numbers, the heat 
dissipation effects tend to becoirie independent of Reynolds number! This 
implies that, contrary to Eq . (~8) , the optimal form of expression for E 

H 

is one which does not contain viscosity v explicitly. This wll] be the 
case in fact if we assume, for example, that 



f \ 

X 


' ■ (^] 


[\] 





where is now hopefully some well behaved empirical function. For 
combining these last two results gives 







(- 10 ) 



which is not explicitly dependent on v. This form is similar to that 
suggested by Prandtl in Pef. (2). 

The fact that the viscous dissipation occurs in a form vhich does not 
explicitly contain the viscosity seems rather paradoxical at first sight, 
and calls for some comment. It can be given a reasonable interpretation by 
considering turbulence energy processes from the spectral point of view. 

It is well known that the work input, term T, occurs mainly at the long 
vjave length end of the spectrum, and that the dissipation to heat, term II, 
occurs mainly at the short wave length end. At high Reynolds numbers, 
there is a broad range of wavelengths in between, the so-called inertial 
range, for vhich both work input and heat dissipation are negligible. Each 
v^ave length within this inertial range simply receives energy from the 
longer wave length components and transmits this energy to the shorter 
wave length com.ponents. It would appear that the rate of transmission of 
energy thru the spectrum is largely controlled by this process in the 
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inertial range, rather than by the viscous dissipation process itself. 
Apparently the viscous dissipation rate easily adjusts itself as may be 
required to dissipate all of the energy coming through the inertial range. 
Since the viscous process is not rate controlling, the effect of viscosity 
tends to drop out of the experimental picture. 

These experimental facts and theoretical considerations suggest 
that the heat dissipation should be analyzed by the same general method 
previously employed in connection with the eddy viscosity. For example, 
we might assume tentatively that 



(” 11 ) 

Then non-dimensionalizing on the basis of a and A as common parameters, we 
obtain 

rAE, 



H 



= p. ^ p. ill Aa V 
^ I A’ v’ A’’ 



(- 12 ) 



which is seen to be equivalent to Eq. (-10) 

This result suggests that is a well behaved universal function 
which can be established empirically. However sound in principle, this 
idea is rather difficult to apply at the present time owing to the paucity 
of suitable experimental data. It does seem plausible that for mean flows 
which are not too strongly unsteady, the effect on 6^ of the parameter 

might prove to he negligible. On the other hand, efforts to establish 



r 



6' as a function of 






alone have not as yet proved conclusive one way 
or the other. This is due in part to limitations in the data and in part 
to limitations in time. 

Instead of following the above line of investigation to the end, we 
discovered the follow^ing more promising modification. Numerical experimenta- 
tion suggested that it might provide a more satisfactory alternative. 
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Recall that the quantities A and i in our ipodel are themselves computed 
from the integrals I and J. Hence it amounts merely to a rearrangement of 
form to rewrite Eq. (-11) as follows: 

(q>ljd,v,y,.*.*) (— 13 ) 

Now, non-dimensional izing on the basis of q and J as common parameters, we 
obtain 






= e 



_(qJ) 



2/3’ 





N 

J 


lj“J ’ 


2 

''a ^ 



1/3 



(-14) 



Actually, the significance of the dimensionless group on the left was 
established by certain somewhat speculative physical arguments v^hich will 
not be repeated here. (See ref. 11) It was ‘found, however, on the basis of 
numerical investigation, that this particular dissipation parameter 6 seems 
to remain more nearly constant, and to agree better with physical data 
than does the parameter resulting from the earlier analysis. Hence 
this is the form that was retained for our present heuristic turbulence 
model. 



Here (?> is a slowly varying dimensionless function. The dependence 

5 ^ 1/3 



of B on the quantities 



(qJ) 



^ and 
2/3 V 



J 



V7as not actually investigated. 



due to time limitations. Hov^ever, it was found heuristically that 3 seems 



to the nearest vail, if any. 



to depend only on the relative distance 
and that far from any walls B becomes essentially constant. 

The following empirical equation seems to produce reasonably good 
curve fits, namely. 



e = 



3.7 



1 e 






(-15) 
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Some data which tends to confirm this formulation is given in 
Appendix C. In any case, it is clear that the energy dissipation to heat 
needs further investigation. 

Next consider term IV of the basic energy equation. It is useful 
for purposes of discussion 



i 



P + 



u .u 
-^,J. 



= U.p + u. , 

1 1 2 



u .u 



(- 16 ) 



The second of these terms clearly represents the net transport of the 

(U ,u^ 

fluctuating turbulent kinetic energy 



by the action of the turbulent 

is 



u .u . 



velocity fluctuations u^. Clearly the turbulent kinetic energy j 
both a scalar and an extensive property. It may be regarded as something 
which is physically transported along with the element of fluid mass. In 
this respect it is quite analogous to other extensive properties, such as 
salinity, for example. It is customary and appropriate to represent such 
net turbulent transport by an appropriate empirical diffusion coefficient 
Hence we write 



ru ,u 

_J_ .3 


- -r ’ 






- -r ' 


faE ' 


2 . 


— — C 


9x. 1 
1 


1 2 . 


~ L 


9x . 



(-17) 



On the other hand the pressure velocity correlation term u^p is in 
a somewhat different category owing to the fact that pressure p is an 
intensive property. It is not an extensive property which can be regarded 
as being transported along with the fluid mass. Me note however that p is 
a dependent variable. This means that whenever the spatial distribution 
of the velocity fluctuations u^ is specified, the corresponding spatial 
distribution of p is likewise fixed. Hence, looking at the matter statisti- 
cally, ve can say that p is correlated in some fashion with u^. In fact, 
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2^1 . For example, 



taking into account dimensional considerations, we can Improve this state- 
ment and say that p is statistically correlated with 

in the limiting case of an inviscid fluid, Bernoulli’s equation informs 
us that regions of higher than average kinetic energy x^^ill tend to be 
regions of lower than average pressure, and vice versa. Consequently we 
are justified in writing the pressure term in the analogous form 



u^p = -e 



3E 



9x . 
1 



(- 18 ) 



Novr the overall effect becomes 



“i('’ -"-H) ■ "''(H:) 



(- 19 ) 



Normally e’ is a positive quantity. Hov^ever, because of the generally 
negative correlation between pressure and kinetic energy, we expect e" to 
be usually negative and smaller than e’. Hence (e*' + e’) should be 
positive but smaller than e’ alone. Moreover, e” and e* are both of the 
same dimension as the kinematic eddy viscosity e. Also e’ can be expected 
to be roughly of the same order of magnitude as e. For these reasons, it 
is advantageous to write. 



(e*' + e’) = ye (-20) 

where y is not a dimensionless coefficient roughly of order of magnitude 
uni ty . 

Hence finally v/e write 




Various numerical experiments suggest that the dimensionless em.pirical 
coefficient like the coefficients a and F, depends only on the distance 
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y to the nearest fixed boundary, if any. 
found to give reasonably good results. 



Y 



1.4 



0 . 4e 




^2 

1 



The following expression has been 



(- 22 ) 



Turning next to term V of the basic energy equation, the molecular 
diffusion effect, it may be said that this will normally be very much 
smaller than the other terms in the equation. It may safely be neglected 
everywhere except in the very thin laminar sublayer region along any solid 
boundary. We shall ignore it in our mode]. 

Finally ve combine the foregoing steps to obtain an overall heuristic 
energy equation as follows. Let 

= 2E (-23) 

Then the energy equation becomes 



8E 



8t 



i i L 



8E 



9x , 

1 ' J 



(-24) 



where e, B, and a are fixed by the following empirical expressions, namely, 




Y = 1.4 - 0.4e 



(-25) 



(-26) 



(-27) 
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/APPENDIX A 



Derivation of Exact Energy Equation for Turbulence 
E quation of r>otion 

Let prime marks on a symbol be used to denote the sum of the mean 
value of that quantity plus the turbulent perturbation. That is, 3et 



Also let 



u : = u . + u . 
j j j 



p’ = p + p 



T : . = T . . + t . . 

ij ij 



denote viscous shear stresses. 

The equations of continuity and motion now become 



3u: 



oX . 
1 



= 0 



(^-1) 



(- 2 ) 



(-3) 





fST! .1 


[st ] ■ 


Ui ] 



8x . 11 

2 



(-4) 



Tl;e corresponding energy eouation is found by multiplying the equation 



of motion by Ul. Thus 
J 



. rU’.US rJT' 

= u* — ^ 

2 ] 



- u: (I'lU!) 

J 9x^ 1 J 



(-5) 



Now we substitute Eqs. (-1) and (-2) into (-3) , (-4) and (-5) , then 
expand, average and simplify. All terms linear in perturbation quantities 
vanish in the averaging. Therefore, the averaged equations of continuity, 
motion and energy become, respectively, 
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(-7) 



(- 8 ) 



Also, multiplying tl>e averaged equation of motion (-7) through by l\ 
gives the energy equation for the mean flov:, namely. 



fV.V A 


r3T. .’I 


LJ_1 


= 0 


1 2 J 


J 



(-9) 



Now subtracting (-1C) from (-9) gives the energy equation for the 
turbulence itself, that is 
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It is convenient to symbolize the mean turbulent kinetic energy, the 
instantaneous kinetic energy, and the Reynolds stresses, respectively, by 



u ,u . 

F = 

2 

u ,u , 

E* = 



(- 11 ) 

(- 12 ) 



T. . = -U.U. 

IJ 1 J 



(-13) 



Therefore the energy equation becomes 
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J ■ "j 


[3x. J 



r3U. 



T — (u.F,') 
3x. 1 



(- 14 ) 
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The viscous term can novj be developed further as follows. Consider 
tlie identity 



u . 



Jt 



ii 



5 /'3u. 



j I 8x . 

Now the viscous stress perturbation may be written 



where 



t . . = 
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‘ ij 'ij 
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Therefore 
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= 0 + V Y . , 



Substituting (-19) and (-20) into (-15) gives 
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at. . 
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Substituting this back into (-14) and rearranging, we find 
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(- 15 ) 
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(-1B) 
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(- 20 ) 



(- 21 ) 



(- 22 ) 
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Two further simplifications are advantageous. Recall that 






3U. 

— J + — 1 
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■ i J 



(- 23 ) 



Noting further that and are both symmetrical tensors, we write 



or 



T. .r. . = T, . 
ij ij 
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By an exactly analogous argument we find that 
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Consequently, the turbulence energy equation becomes finally. 
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(-27) 



It is important to notice that this equation is complete and exact. 
Each term has a definite physical significance as discussed in the main 
body of the report. The method of evaluating certain of the terms heuris- 
tically is also discussed in the main report. 
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APPENDIX B 



Local Length Scale for Axi-Syimnetric Mean Flows 



Consider the following sequence of relations pertaining to the 
heuristic model, where the various integrals are summed over the entire flow 
field. 



W 

H 



w 



Ax* Ax 



= e = Weighting function 

= /Wdv’ = Normalizing integral 
W 



H 



= Normalized weighting function 



(B-1) 

(- 2 ) 

(-3) 



I 

J 

X 



2 

2 

2 



= /i^fi^dv’ = First characteristic integral 

2 

= /w(r^Q^) dv’ - Second characteristic integral 



= Local scale parameter 



(-4) 



(-5) 

(- 6 ) 



For axi“S>Tnmetric flov:s, the volume integrals (-2), (-4), and (-5) 
reduce to surface integrals over the ireridlonal plane. 

In this case the volume element may he written in the form 
dv’ = r’de^drMz’ 

Also in cartesian coordinates 

Ax-Ax = (x^ - x)^ + (y’ - y)"" + (z* - z)^ (-7) 

We may express the coordinates of the fixed point x and of the 
variable point x’ in terms of cylindrical coordinates as follows. 
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X = 0 



x' = r' sin 0' 



y = r y ' = r ' cos 6 ' (-8) 

z = z z ' = z ' 

Upon substituting (-8) into (-7) and rearranging, we can write 

Ax*Ax = + 2rr'(l - cos 0') (-9) 

where 

^ = (r’ - r)^ + (z' - z)^ (-10) 

2 

Mow let us consider the typical integral I as defined by Eq . (-^) . 

By making use of the foregoing relations, this integral can be reduced to 
the required form. In addition, it is convenient to utilize the following 
auxiliary variable, namely, 



C 



2 




(- 11 ) 



Eq. (-4) is now expressible in the form 




2ir 



k I ~ ®'^de'j yr^Ojdr'dz' (-12) 



Now since Q is a function of z’ and r* only, but independent of 0’, 
the above integration with respect to G* can be made once and for all. 
This leads naturally to the definition of the auxiliary function 



^U) 



2-n 



/If 



e O - cos e')^g, 



(-13) 



O 

The function ip defined above must be computed numerically. The 
factor was inserted on the right so as to achieve the simplification 
that at large values of C the function i' approaches the limit unity, that 
is, that 
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JiiTii ij; = 1 

^-KX3 



(- 14 ) 



It then happens that for small values of C 

£im = 2/T c (-15) 

C->o 



The function ijj is shovm in Fig. R-1 . Computed values are listed in 
Tabic B-1. 

2 

Once the function if is knov7n, the integral I may be rewritten as 
follows . 



I 



2 




n^dr ' dz' 



(-16) 



Similarly 









yrr' 



7r’ (n^fi'^)dr'dz' 



(-17) 



The scale parameter X is now found from (-6) . Notice that in 
computing X^ from (-6) the factors ^ (-16) and (-17) cancel from the 

result . 



Eqs. (-16) and (-17) illustrate how the characteristic integrals 
reduce from volume to surface integrals for the axi-symmetric case. Eos. 
(-11) and (-13) define il' which is a function characteristic of axi-symmetric 
mean flows. 
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1.0 (asymptote) 
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FIGURE B-l . GRAPH OF THE FUNCTION W ) 



TABLE B-1 



COMPUTED VALUES OF THE FUNCTION 
(Refer to Eqs. B-11, B-13) 



0.0 


O.OOOC 


0.1 


0.3475 


0.2 


0.6555 


0.3 


0.8955 


0.4 


1.0562 


0.5 


1.1433 


0.6 


1.1739 


0.7 


1.1687 


0.8 


1.1461 


0.9 


1.1187 


1.0 


1.0936 


1.2 


1.0581 


1.4 


1.0387 


1.6 


1.0279 


1.8 


1.0213 


2.0 


1.0169 


2.5 


1.0105 


3.0 


1.0072 


4.0 


1.0040 


5.0 


1.0025 


oo 


1.0000 
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APPENDIX C 



Experimental Data 

Various useful reference items of experimental data are summarized in 
Tables C-1 through C-5 and in Figs. C-1 and C-2. These tabulations and 
diagrams are largely self-explanatory. 

Table C-2 is of special interest in that it shows a comparison 
between theoretical values of dissipation, as computed from the heuristic 
model, and the experimental measurements of Laufer, Ref. (10). The agree- 
ment is satisfactory. 

Laufer 's measured values are knox^m to be too low and he attempted to 
extimate better values as shown. 

Shearing rates were estimated on the basis of a theoretical mean 
velocity profile using Prandtl’s mixing length theory and Nikuradse's 
empirical curve for the mixing length £ , namely, 

I = 0.14 - 0.08r^ - 0.06r^ (C-1) 

The information in Table C-2 is also displayed graphically in Figs. 

C-1 and C-2. There is no particular significance to the fact that the 
ordinate in Fig C-2 is rather than 3 Itself; this merely reflects the 
circumstance that a slightly different notation was being employed at the 
time these calculations were actually made than that which is used in this 
report. 

In view of the uncertainties in the experimental data Itself, the 
agreement between theory and experiment as revealed In Figs. C-1 and C-2 
is regarded as satisfactory. 
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The remaining tables require no particular comment except, perhaps, 
to note that the available data is meagre and not very satisfactory. Note 
that the data given is mostly one dimensional, with slight information in 
the second dimension and none at all in the tliird. There is clearly a 
need for more experimental work. 
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TABLE C-1 



TURBL-LECT ENERGY DATA FOP. A PIPE 
Measurements by Laufer, Ref. (10) 



All quantities non-dlmensionallzed using radius a = 1 and friction 
velocity V* = 1. 



R = 500,000 
e 



r 


2E 




0 


1.65 


2.1 


0.1 


1.80 


2.2 


0.2 


2.05 


2.4 


0.3 


2.45 


3.1 


0.4 


2.90 


3.8 


0.5 


3.50 


4.8 


0.6 


4.20 


6.2 


0.7 


5.00 


8.7 


0.8 


6.00 


13.0 


0.9 


7.40 


22.8 


0.98 


— 






R - 50,000 
e 



2E 




Vii^: 


1.69 


2.15 


— 


1.80 


2.26 


— 


2.06 


2.72 


— 


2.41 


3.28 


1.5 


2.85 


3.96 


2.5 


3.48 


4.76 


3.7 


4.01 


6.23 


5.3 


4.49 


8.10 


7.6 


5.13 


11.3 


11.7 


5.66 


20.9 


22.0 


7.85(peak) 


— 


— 
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TABLE C-2 



MEASURED AlvID COMPUTED ENERGY DISSIPATION 



RATES IN A PIPE AT R = 500,000 

e ’ 



Measurements by Laufer, Ref. (10) 



A1 quantities non-dimensionalized using radius a = 1 and friction 
velocity v* = 1. 



r 






8-^ 










Table C-1 


Calc . 


Eci. (1.3-18) 


Eq. (2.8-14) 


Meas . 


Est . 


0 


1.78 


2.94 


3.70 


1.43 


1 .0 


2.0 


0.1 


2.00 


2.98 


3.70 


1 .61 






0.2 


2.45 


3.09 


4.58 


1.65 


1.2 


2.5 


0.3 


2.95 


3.29 


6.71 


1.45 






0.4 


3.45 


3.62 


7.35 


1.71 


2.0 


4.0 


0.5 


4.30 


4.16 


7.40 


2.42 






0.6 


5.30 


5.03 


7.40 


3.61 


3.0 


6.0 


0.7 


6.70 


6.52 


7.40 


5.91 






0.8 


8.20 


9.58 


7.40 


10.60 


7.0 


13.5 


0.9 


10.00 


18.80 


7.40 


25.40 


15.0 


25.0 


Notes : 


1. Laufer 


regards his 


measured dissipation data 


as too low 


and 



attempts to estimate better balues as indicated. 



1/3 

2. The quantity was obtained from the theoretical 

velocity profile based on the Prandtl-Nikuradse mixing length 
theory using Eg. (C-1) Refer also to Eq. (1.3-7). 



Ill 



4 



S. 

,5Q 

o 

Si 

a> 

> 

c 




0 



□ Loufers Measured Values 

® Laufer's Adjusted Estimates 
Heuristic Model 



0.2 



04 



0.6 



Relative Radius L 
R 



0.8 



1.0 



FIGURE C- 1 . INVERSE DISSIPATION FUNCTION 

IN A PIPE 
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TABLE C-3 



AXIAL VELOCITIES IN A COLD TURBULENT JET 
AT AXIAL STATION ^ = 20 



Measurement by Corrsln and Uberol, Ref. (7) 



V = mean axial velocity 

V = axial velocity at axis 
o 

V = root mean square axial perturbation 

V 1 

bj <2 “ radius to point where — = 



z = axial coordinate 

D = diameter of jet 




0 

0.2 

0.4 

0.6 

0.8 

1.0 

1.2 

1.4 
1.6 
1.8 
2.0 
2.2 

2.4 

2.6 







fv ] 


V 




'• O'* 



1.0 


0.205 


0.95 


0.21 


0.85 


0.215 


0.74 


0.215 


0.625 


0.20 


0.500 


0.18 


0.40 


0.165 


0.32 


0.145 


0.25 


0.12 


0.18 


0.09 


0.135 


0.07 


0.08 


0.045 


0.055 


0.03 


0.04 


0.02 
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TABLE C-A 



AXIAL AND PADIAL VELOCITIES IN A HEATED TURBULENT JET 
AT AXIAL STATION •§ " 15 

Measurements by Corrsln and Uberoi, Ref. (7) 



u = root mean square radial perturbation 
Other notation as in Table C-3 




0 


1.0 


0.17 


0.175 


0.2 


0.97 


0.17 


0.17 


O.A 


0.88 


0.17 


0.16 


0.6 


0.76 


0.17 


0.14 


0.8 


0.62 


0.17 


0.12 


1.0 


0.50 


0.15 


0.09 


1.2 


0.38 


0.13 


0.075 


l.A 


0.28 


0.10 


0.05 


1.6 


0.205 


0.075 


0.035 


1.8 


0.145 


0.05 


0.022 


2.0 


0.10 


0.03 


0.013 


2.2 


0.055 


0.02 


0.005 


2.4 


0.025 


0.01 


0 
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PS u 



TABLE C-5 



AXIAL VELOCITIES IN A TURBULENT JET 
AT TWO AXIAL STATIONS 
Measurements by Laurence, Pef. (8) 
Notation as in Tables C-3 and C-4 



= n 



0 



0.5 



1.0 



1.5 



2.0 



2.5 

3.0 



3.5 



4.0 

4.5 



5.0 



6.0 



7.0 



8.0 



At 




16 





0.43 



0.40 



0.31 



0.22 



0.15 



0.07 




0.082 

0.082 

0.082 

0.084 

0.089 

0.091 

0.090 

0.084 

0.074 

0.061 

0.047 



V 

0.36 

0.34 

0.30 

0.25 

0.19 

0.135 

0.08 

0.025 

0 



■V 

V 



0.073 

0.073 

0.073 

0.074 

0.075 

0.075 

0.075 

0.075 

0.075 

0.072 

0.067 
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